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ABSTRACT 

Context. Beryllium is a pure product of cosmic ray spallation. This implies a relatively simple evolution in time of the beryllium 
abundance and suggests its use as a time-like observable. 

Aims. Our goal is to derive abundances of Be in a sample of 90 stars, the largest sample of halo and thick disk stars analyzed to date. 
We study the evolution of Be in the early Galaxy and its dependence on kinematic and orbital parameters, and investigate its use as a 
cosmochronometer Abundances of Be, Fe, and a-elements of 73 stars are employed to study the formation of the halo and the thick 
disk of the Galaxy. 

Methods. Beryllium abundances are determined from high-resolution, high signal-to-noise UVES spectra with spectrum synthesis. 
Atmospheric parameters and abundances of a-elements are adopted from the literature. Lithium abundances are used to eliminate 
mixed stars from the sample. The properties of halo and thick disk stars are investigated in diagrams of log(Be/H) vs. [ff/H], log(Be/H) 
vs. [Fe/H], and [a/Fe] vs. log(Be/H) and with orbital and kinematic parameters. 

Results. We present our observational results in various diagrams. (0 In a log(Be/H) vs. [Fe/H] diagram we find a marginal statistical 
detection of a real scatter, above what is expected from measurement errors, with a larger scatter among halo stars. The detection of 
the scatter is further supported by the existence of pairs of stars with identical atmospheric parameters and dilferent Be abundances. 
(ii) In a log(Be/H) vs. [a/Fe] diagram, the halo stars separate into two components; one is consistent with predictions of evolutionary 
models, while the other has too high a and Be abundances and is chemically indistinguishable from thick disk stars. This suggests that 
the halo is not a single uniform population where a clear age-metallicity relation can be defined. (Hi) In diagrams of Ri^j,, vs. [ff/Fe] 
and log(Be/H), the thick disk stars show a possible decrease in [ff/Fe] with Rmin, whereas no dependence of Be with R„,in is seen. This 
anticorrelation suggests that the star formation rate was lower in the outer regions of the thick disk, pointing towards an inside-out 
formation. The lack of correlation for Be indicates that it is insensitive to the local conditions of star formation. 

Key words, stars: abundances - stars: late-type - Galaxy: halo - Galaxy: thick disk 



1 . Introduction dependence of the Be abundance with metallicity. However, ob- 

servations of Be in metal-poor stars (Rebolo et al. 1988 ; Gilmore 

. ^ The nucleosynthetic origin of beryllium is different from most ^gg^ ^^j^^.^ ^1 [19973 Boesgaard et al. 1999) find a 

chemical elements. Beryllium is neither a product of stellar nu- ^ j ^j^^^ ^^^^^^^ log(Be/H) and [Fe/H], and 

cleosynthesis nor expected to be created by the standard ho- -^^^ j^jj j^^^ log(Be/H) and [O/Hfl Such slopes ar- 

mogeneous pnmordial nucleosynthesis in a detectable amount ^j^^^ ^^^^^^^ ^ ■ ^^^^^^^ .^^ production 

(Thomas et al. 1993 ). Its single long-lived isotope, ^Be, is a pure mechanism is independent of ISM metallicity. Thus, the dom- 

product of cosmic-ray spallation of heavyO^iosdy CNO) nuclei -^^^^ production mechanism is now thought to be the collision 

indie interstellar medium (Reeves et al. [MOl Meneguzzi et al. cosmic-rays composed of accelerated CNO nuclei with pro- 

ilSjJi- tons and a-particles of the ISM (Duncan et al. |1992t Casse et al. 

Early theoretical models of Be production in the Galaxy as- \ igg^ [ Vangioni-Flam et al | 1998| l 
sumed the cosmic-ray composition to be similar to the compo- 

sition of the interstellar medium (ISM). In this scenario. Be is ^n the other hand, in hght of some controversyabout the be- 
produced by accelerated protons and a-particles colliding with Savior of O in metal-poor stars. Fields & Olive ([19parguethat 
CNO nuclei of the ISM (Meneguzzi & Reeves |T975l Vangioni- ^ secondary behavior cannot be excluded. King dMI] |2Q02j, 
Flam et al. 119901 Prantzos et al. |1993, resulting in a quadratic however, reassessing data of different oxygen indicators, as well 
as data for other a-elements, show that a metal-poor primary 



* Based on observations made with ESO VLT at Paranal mechanism is necessary to explain the observations. 
Observatory, under programs 076.B-0133 and 077.B-0507, and on data 

obtained from the ESO/ST-ECF Science Archive Facility and the UVES 

Paranal Observatory Project 266.D-5655. ' [A/B] = log [N(A)/N(B)]j, - log [N(A)/N(B)]o 
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As a primary element and considering cosmic-rays to be 
globally transported across the Galaxy, one may expect the Be 
abundance to be rather homogeneous at a given time in the early 
Galaxy. It should have a smaller scatter than the products of stel- 
lar nucleosynthesis (Suzuki et al. 1999; Suzuki & Yoshii EOOTb . 
Thus, Be would show a good correlation with time and could 
be employed as a cosmochronometer for the early stages of the 
Galaxy (Beers et al 2000, Suzu ki & Yoshii 2001} . 

Pasquini et al. (120041 120071 1 have tested this suggestion de- 
riving Be abundances in turn-off stars in the globular clusters 
NGC 6397 and NGC 6752. The Be ages de rived f rom a model 
of the evolution of Be with time ( Valle et al. I2002I I are in excel- 
lent agreement with those derived from theoretical isochrones. 
Moreover, the Be abundances of these globular cluster stars are 
similar to the abundances of field stars with the same metallicity. 
These results strongly suggest that the stellar Be abundances are 
independent of the environment where the star was formed and 
support the use of Be as a cosmochronometer. 

Using abundances determined by Boesgaard et al. d 19991 1. 
Pasquini et al. (I2005I I extended the use of Be as a time scale 
to a sample of 20 halo and thick disk stars and investigated the 
evolution of the star formation rate in the early-Galaxy. Stars be- 
longing to the two different kinematic components identified by 
Gratton et al. (I2003al l seem to separate in a log(Be/H) vs. [O/Fe] 
diagram. Such separation is interpreted as indicating the forma- 
tion of the two components to occur under different conditions 
and time scales. 

Our aim in this work is to better understand the evolution 
of Be in the early Galaxy and its dependence on different pa- 
rameters, in particular on the stellar population. We analyze an 
unprecedentedly large sample of halo and thick disk stars, and 
further investigate the use of Be as a cosmochronometer and 
its role as a discriminator of different stellar populations in the 
Galaxy. The sample and the observational data are described in 
Sect. |2] The details on the adopted atmospheric parameters are 
given in Sect. |3] The abundance determination and its uncertain- 
ties are discussed in Sect. |4] while a comparison with previous 
results of the literature is presented in Sect. |5] Lithium abun- 
dances are used to clean the sample from mixed stars in Sect. |6] 
The Galactic evolution of Be is discussed in Sect. [T] the use of 
Be as a chronometer is discussed in Sect. [8] while a summary is 
given in Sect. |9] 

2. Sample selection and observational data 

The sample stars were selected from the compilation by Venn 
et al. (I2004I I of several abundance and kinematic analyses of 
Galactic stars available in the literature. Using the available kine- 
matic data, Venn et al. (,2004) calculated the probabilities that 
each star belongs to the thin disk, the thick disk, or the halo. A to- 
tal of 90 stars were selected for this work: 9 of them have higher 
probability of being thin disk stars, 30 of being thick disk stars, 
49 of being halo stars; and 2 have 50% probabilities of being halo 
or thick disk stars. We simply assume a star belongs to a certain 
kinematic group when the probability of belonging to that group 
is higher than the probability of being in the other two groups. 
One of our aims is to compare stars of different populations but 
of similar abundances, so we tried to maximize the metallicity 
overlap between the two sub-samples. The halo stars range from 
[Fe/H] = -2.48 to -0.50 and the thick disk stars from [Fe/H] 
= -1.70 to -0.50, although strongly concentrated in [Fe/H] > 
-1.00. The sample stars are listed in Table [1] together with V 
band magnitudes, parallaxes, absolute magnitudes, bolometric 
corrections, bolometric magnitudes, luminosities, and informa- 



tion on multiplicity. Details on these are given in the following 
sections. 

Spectra for all stars were obtained using UVES, the 
Ultraviolet and Visual Echelle Spectrograph (Dekker et al. 200^ 
fed by UT2 of the VLT. UVES is a cross -dispersed echelle spec- 
trograph able to obtain spectra from the atmospheric cut-off at 
300nmto ~ llOOnm. 

For 55 stars, new spectra were obtained in service mode dur- 
ing two observing periods between October 2005 and September 
2006. UVES was operated in dichroic mode with cross dis- 
persers #1, #3, and #4 resulting in spectral coverage of AA 300- 
390 nm in the blue arm and AA 420-680 nm and AA 660-1100 
nm in the red arm, with some gaps. The spectra have a resolv- 
ing power of R ~ 35 000 and S/N between 25 and 150 in the Be 
region. The reduction was conducted with the UVES pipeline 
within MIDAS. While the blue arm of the spectra was always 
found to have enough quality for the analysis, the same was 
not true for the red arm that, in some cases, presented resid- 
ual fringing. For them, a new reduction of the raw frames was 
conducted using the most recent release of the UVES Common 
Pipeline Library (CPL) recipes within ESORex, the ESO Recipe 
Execution Tool. 

Reduced spectra for 4 stars were downloaded from the 
UVES Paranal Observatory Project (POP) Library - a library of 
high-resolution spectra of stars across the HR diagram (Bagnulo 
et al. 2003 ). In this library the stars were observed with both 
dichroic #1 and #2 covering almost all the interval between 300 
and 1000 nm. The spectra have R ~ 80000 and a typical S/N 
ratio varying from 60 to 400 in the Be region. 

For the remaining 3 1 stars, archive raw data in the AA 300- 
390 nm region were retrieved from the ESO/ST-ECF science 
archive facility and employed in the analysis. For some stars, 
in particular those that do not have previous determinations of 
lithium abundance in the literature, raw data in the AA 420-680 
nm region were also retrieved. The spectra were reduced using 
the UVES pipeline within MIDAS and ESORex. The resolving 
power of these spectra varies between 35 000 and 50 000 and the 
S/N ratio varies between 45 and 170 in the Be region. The log 
book of the observations is given in Table |2]in the Appendix. 

3. Atmospheric parameters 

3.1. Source of the parameters 

All the selected stars have been targets of previous abundance 
analyses. In this work, we decided to adopt the atmospheric pa- 
rameters, effective temperature (Tes), surface gravity (log g), mi- 
croturbulence velocity (^), and metallicity ([Fe/H]), determined 
by previous works. The adopted parameters are given in Table[3] 
We adopted as main reference the same high-resolution analy- 
ses quoted in the compilation by Venn et al. (1200 4'). Our sample 
contains 59 stars analyzed by Fulbright ( 120001 FOO hereafter), 15 
by Nissen & Schuster ( 1997 NS97 hereafter), 7 by Edvardsson 
et al. ( 1993, Ed93 hereafter), 4 by Stephe ns & B oesgaard ( 20021 
SB02 hereafter), 2 by Prochaska et al. (|2000l PrOO hereafter), 
1 by Gratt on & S neden ( [19881 GS88 hereafter), 1 by Gratton 
& Sneden fTWTl GS91 hereafter), and 1 by Mc William et al. 
( 119951 McW95 hereafter). 

Even though different methods were adopted to determine 
the parameters, we do not expect this choice to introduce any 
systematic effect in our analysis. The likely effect of adopting 
atmospheric parameters with possibly different scales is an in- 
crease in the scatter of the abundances. We note that Venn et al. 
(120041 ) argue that a superficial analysis of their whole sample did 
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Reference 


Temp. 


Surface gravity 


Microt. 


Model atm. 


LTE/NLTE? 


FOO 


Exc. Eq. of Fe i lines 


Ion. Eq. Fe i and Fe ii 


Fe I lines 


Kurucz + overshooting 


LTE 




log (EW/A) < -4.80 


log (EWM) < -4.80 








SB02 


Exc. Eq. of Fe i lines 


Average of Ion. Eq. of 


Fe I lines 


Kurucz 


LTE 




log (EWM)<-5.15 


Fe i/Fe ii and Ti i/Ti ii 


agrees with Ed93 relation 


Castelli et al. ('1997') 




NS97 


Exc. Eq. of Fei lines 


Ion. Eq. Fe i and Fe ii 


Fe I lines 


(NEW)MARCS (Ed93) 


LTE 




log (EW//1) < -5.15 




with;^' > 4 eV 






Ed93 


Stromgren photometry 


Balmer discontinuity 


Empiric relation 


(NEW)MARCS 


LTE 




calib. of (b - y) and /J 


index, Ci 


dependent on T^g and log g 







not show a major inconsistency or zero point difference between 
the results. We refer the reader to the original papers for detailed 
descriptions of the methods and comparisons between the results 
with other literature determinations. We note, however, that all 
the papers claim reasonable agreement with previous determina- 
tions. Moreover, as the parameters of most of the stars analyzed 
here (59 out of 90) come from a single reference, FOO, we ex- 
pect any possible difference in the scale of the parameters to 
have only a minor effect on the abundances. 

For the completeness of the discussion, we list some infor- 
mation in Table|?]on the methods adopted by each relevant work 
as the reference for the atmospheric parameters. As is clear from 
TableH FOO, NS97, and SB02 adopt very similar methods based 
on spectroscopy to determine the atmospheric parameters. On 
the other hand, the method adopted by Ed93 is completely di- 
verse, based only on photometry. In NS97 the same model at- 
mospheres used by Edvardsson et al. are adopted. In addition, 
their work is a differential one with respect to two bright stars 
analyzed by Ed93. They therefore obtain parameters essentially 
on the same scale as the ones derived by Ed93, as demonstrated 
by their Fig. 03, which shows excellent agreement between their 
Teff with temperatures derived using (b - y) calibrations. 

As Be abundances are calculated from lines of the ionized 
species, log g is the most relevant parameter for the analysis. 
Our sample stars are relatively bright and nearby, so we used 
Hipparcos parallaxes (lESA 1997 ) to estimate gravities. Apart 
from eight stars that show a significant difference (larger than 
0.28 dex), the agreement is excellent. A linear fit to the points 
yields a line that is statistically indistinguishable from the x = y 
line, with a r.m.s. of 0.20 dex when all the points are considered, 
or of 0.14 dex when the discrepant stars are excluded. These 
eight stars have uncertain parallaxes o-Jtt > 0.20 or are binary 
stars. 

Thus, only a comparison between the results of FOO and 
Ed93 is necessary to ensure the consistency between the atmo- 
spheric parameters of most of the sample stars. Such a compari- 
son is shown below. 

3.2. Comparison between Edvardsson et al. and Fulbright 

The comparison between the results of FOO and Ed93 was con- 
ducted to assure the detailed understanding and proper identifi- 
cation of any possible systematic effect introduced in the analy- 
sis. 



^ We do not include in Table |4] details on the analyses by Gratton & 
Sneden ( 1988 , 1991 ) and by McWilliam et al. ( 1995 ), which contribute 
only 1 star each, and on the one by Prochaska et al. ( 2000) , with only 2 
stars, since they will have no influence at all on a possibly larger scatter 
of the results. 



Table 5. Uncertainties in the adopted atmospheric parameters. 



Ref. 


Teff (K) 


logg 


f (kins"') 


[Fe/H] 


FOO 


±40 


±0.06 


±0.11 


±0.04 


SB02 


±75 


±0.35 


±0.30 


±0.10 


NS97 


±50 


±0.10 




±0.05 


Ed93 


+ 100 

-50 


±0.20 


±0.30 


±0.10 


adopted 


±100 


±0.15 


±0.30 


±0.15 



Ten stars are shared between the two analyses. Of these, 
seven are included in our sample. The atmospheric parameters 
of all the ten stars derived by the two papers are compared in 
Fig. [T] The effective temperatures show excellent agreement, 
with a mean difference of 46 K and a maximum difference of 
85 K. Although the values derived by FOO are systematically 
lower than the ones derived by Ed93, the difference is small, 
well within the uncertainties of the calculations. The difference 
in surface gravity is also small, with a mean of 0.11 dex and 
a maximum of 0.19 dex, and might be related to the different 
adopted temperature scales. The differences in gravity are close 
to the uncertainty value. Finally, the metallicity ([Fe/H]) values 
show a mean difference of 0.10 dex and maximum of 0.21 dex. 
This difference seems to be systematic, because the values of 
FOO are always slightly higher than the ones derived by Ed93. 

The comparison shows that the differences between FOO and 
Ed93 are not large, but can be ascribed mostly to the uncertain- 
ties in each analysis, except for the microturbulence (Fig. [TJ. 
This shows that no strong systematic effect will be present on 
our abundances and argues in favor of our assumption that the 
scales of the different papers are indeed similar, and the use of 
the combined set of parameters introduces only a small and ac- 
ceptable dispersion on the abundances. However, we note that 
the observed difference in [Fe/H] might introduce some disper- 
sion on the data points when analyzing plots like log(Be/H) vs. 
[Fe/H] or [Be/Fe] vs. [Fe/H]. This point has to be kept in mind 
and will be recalled when discussing our results. 

3.3. Uncertainties of tiie atmospiieric parameters 

In Table|5]we list the uncertainties in the atmospheric parameters 
as calculated by each of the four papers that contribute more than 
two stars to our sample. 

As seen in Table |5] the adopted works claim a rather low 
unceitainty for their Teff values, ranging from 40 to 100 K. As 
discussed before, for the ten stars in common, the mean differ- 
ence between the temperatures derived by Ed93 and FOO is 46 
K with maximum of 85 K. These values argue in favor of the 
quoted uncertainties. Nevertheless, in this case we decided to 
adopt a rather conservative value of 100 K as representative of 
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Fig. 1. Comparison of the adopted atmospheric parameters for stars shared by FOO and ED93. The error bars with the values Hsted 
in the original papers are also shown. 



the associated uncertainty. This value equals the higher of the 
values suggested by Ed93 and might be a good estimate of the 
effect of mixing the temperature scales of the different papers 
together. However, as we will see in the following sections the 
uncertainty in Teff has little influence on the derived beryllium 
abundance. 

For the log g values, the adopted papers claim uncertainties 
ranging form 0.06 to 0.35 dex. The comparison between Ed93 
and FOO shows a mean difference of 0.11 dex with maximum 
of 0.19 dex. Moreover, the comparison with log g values esti- 
mated using parallax show an r.m.s. of 0.14 dex. In particular, 
this comparison seems to show that the log g values are rather 
well constrained, in spite of the larger uncertainties quoted by 
Ed93 and SB02. Therefore, in this case we adopted a value of 
0.15 dex as representative of the uncertainty in log g. 

The quoted uncertainty on the microturbulence velocity (^) 
varies from 0.1 1 to 0.30 km s"' in the adopted papers. The com- 
parison between Ed93 and FOO, however, shows larger differ- 
ences with a mean of 0.31 km s ' and maximum of 0.54 km 
s~^ In this case, however, we simply decided to adopt the higher 
value of the quoted uncertainties, 0.30 km s"' , as representative. 
Again, as we see below, this choice and the large differences in 
^ have no significant impact on the beryllium abundances. 



The adopted papers list values for the metallicity ([Fe/H]) 
ranging from 0.04 to 0.10 dex for its uncertainty. The compari- 
son between Ed93 and FOO, however, shows a mean difference 
of 0.10 dex with maximum of 0.21 dex. Mainly because of this 
comparison, we decided to adopt a somewhat higher value as 
representative of the uncertainty than the ones quoted in the 
adopted papers, 0.15 dex. As for the other parameters, except 
for log g, its influence on the final beryllium abundances is weak. 
The adopted uncertainty values are listed in Table |5] 



4. Abundances determination 

4.1. Synthetic spectra 

Abundances of beryllium and, for some stars, of lithium were 
calculated using synthetic spectra. In this work, we adopted the 
codes for calculating synthetic spectrum described by Barbuy 
et al. ([2003] l and Coelho et al. ("2005). We use the grids of 
model atmospheres calculated by the ATLAS9 program (Castelli 
& Kurucz '2003 ), without overshooting. These models assume 
local thermodynamic equilibrium, plane-parallel geometry, and 
hydrostatic equilibrium. Codes for interpolating among the grids 
were adopted. 
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3129.5 3130 3130.5 3131 3131.5 



Fig. 2. Fit to the region of the beryllium lines in the solar spec- 
trum and in the stars HIP 105888 and HIP 24030. The crosses 
represent the observed spectrum, the solid line the best synthetic 
fit, and the dashed lines represent changes of +0.20 dex in the 
Be abundance of the best fit. A solar abundance of A(Be) = 1.10 
was determined. 

On the spectra of late-type stars, only the Ben - ^Pq res- 
onance lines at 3131.065 A and 3130.420 A are strong enough 
to be useful for an abundance analysis. These lines are near the 
atmospheric cut-off at 3000 A in a region of low detector sensi- 
tivity. This near-UV region is extremely crowded, full of atomic 
and molecular lines, some of them still lacking proper identifica- 
tion. The determination of Be abundances, thus, needs to be done 
with spectrum synthesis taking all the blending nearby features 
into account. Our database of molecular lines include the OH 
(A^E-X^n), NH (A^n-X^S), and CN (B ^H-X^S) band systems 
as implemented by CastiUio et al. d 19991 1 and the CH band sys- 
tems (A^A-X^n), ( B^E-X ^n), and (C^S-X^H) as implemented 
by Melendez et al. ( I2003I I. all of them affecting this spectral re- 
gion. 

4.2. Beryllium 

The line list of atomic lines compiled by Primas et al. d 19971 1 
was adopted in this work. This same line list was adopted in 
many analyses of Be abundances in literature (Garcia Perez & 
Primas 2006 Primas et al. T9971 12000ai l2000b Randich et al. 
12002 2007). The adopted log gf of the Be lines are those com- 
monly used in the literature, -0.168 and -0.468 for 3131.065 
A and 3130.420 A, respectively. In this line list, the blending 
line affecting the blue wing of the Be 3131 line is assumed to be 
an Fei line in 3131.043 A, with log gf = -2.517 and x = 2.85 
eV. The parameters of this line were constrained using several 




3129.5 3130 3130.5 3131 3131.5 

m 



Fig. 3. Fit to the region of the beryllium lines in the stars HIP 
24316, HIP 1 15167, and HIP 48152. The lines are as in Fig.|2] 



stars of different parameters and metallicities (see Primas et al. 
IT997 for details). The proper identification of this line, however, 
remains controversial in the literature. Castilho et al. (| 199 91), for 
example, adopt an Fe i line at 3 1 30.995 A in the analysis of Be in 
red giants, with log gf - -3.30 wdx - 3.00 eV. Other possibil- 
ities explored in the literature include an Mm line at 3131.037 
A adopted by Garcia Lopez et al. (1995a) and an Mnii line at 
3131.015 adopted by Boesgaard & King ( [2002ll. In this work, 
we opted for the Fe i included by Primas et al. ( 119971 ), in partic- 
ular to maintain the consistency with all the papers that adopted 
the same atomic line list. 

There is currently a controversy about a possible missing 
continuum opacity source affecting the near-UV region, hence 
Be abundances (Asplund 120041 and references therein). The 
code used to calculate the synthetic spectra does not include 
bound-free opacities due to metals (Smiljanic & Barbuy 2007ll. 
However, we are mostly interested in relative abundances be- 
tween stars of similar composition, so the bulk of the analysis is 
not affected by this uncertainty. The slope between Be and Fe or 
a could be affected. With respect to this, we note that our most 
metal-rich stars have [Fe/H] - -0.50; for metal-poor objects, the 
missing opacity effect is expected to be negligible. 

In fact, Balachandran & Bell (1998) argue that the Fei 
bound-free opacity should be increased by a factor 1.6. This is 
equivalent to an increase of 0.20 dex in Fe abundance. We ran a 
test using the Kurucz suite of codes to test the influence of this 
increase in Fe abundance in the calculation of Be abundance for 
a model with [Fe/H] = -0.50, appropriate for star HIP 31639. 
The difference of 0.022 dex is negligible. We therefore expect 
that this issue does not really affect our analysis and conclusions 
in any significant way. 
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With the line list described above, we fitted the Be lines re- 
gion in the solar UVESQ spectrum, adopting the parameters Teff 
= 5777 K, log g = 4.44, and ^ = 1 .00 km/s, and obtained A(Be) = 
1.10 (Fig.|2|. This abundance is in excellent agreement with the 
one found by Chmielewski et al. ( |19751 l, A(Be) = 1.15, usually 
adopted as the reference photospheric solar abundance. We also 
note excellent agreement with the abundance derived by Randich 
et al. (I2002I I. A(Be) - 1.11, using this same atomic line list but a 
different molecular line list and by fitting the Kurucz Solar Flux 
Atlas (Kurucz et al. |1984| l. 

From the 90 stars in the sample, abundances of Be were de- 
termined for 83: 76 detections and 7 upper limits. The 7 stars 
for which abundances could not be calculated were affected by 
different problems, some with profiles affected by cosmic rays, 
and others were found to be spectroscopic binaries. These stars 
will not be considered further in the analysis. Examples of the 
fits used to derive the Be abundance are given in Figs. |2]and[3] 

For most of the stars, both beryllium lines are well-fitted with 
the same abundance. In some stars, however, the lines are best 
fitted with slightly different values. In most cases, the difference 
amounts to ~ 0.04 dex, a value that can be completely attributed 
to the uncertainty of the fitting procedure itself and/or to the de- 
termination of the continuum. In one extreme case, however, the 
difference amounts to 0.24 dex (HIP 31639). This discrepancy 
is likely due to an inefficient treatment of the blending features 
caused, for example, by abundance ratios that are significantly 
different from solar of the blending lines directly affecting the 
Be lines. For these discrepant cases, we consider the abundance 
derived solely from the 3131 A line to be more reliable, since its 
profile is visually less affected by blends. In Table[3]we list both 
the beryllium abundance determined solely from the 3 1 3 1 A line 
and the average of the two lines. 

We tested the dependence on atmospheric models by also 
using a version of the Kurucz model atmospheres with enhanced 
abundances of a-elements, [c/Fe] = -1-0.40 dex. Metal-poor stars 
are known to have enhanced abundances of a elements. Some of 
these elements are particularly important as electron donors and, 
as such, may influence some continuum opacity sources (bound- 
free and free-free absorption of H , for example). We found, 
however, this effect to be small. In the following discussion, we 
adopt the beryllium abundance derived solely from the 3131 line 
with normal model atmosphere as the reference value. 

4.3. Lithium 

In our analysis, we expect to recover the initial beryllium abun- 
dance the star had at the time of its formation, since this is the 
one produced by the interaction between the Galactic cosmic- 
rays (GCR) and the ISM. It is important then to identify stars 
where significant mixing of the photospheric material with deep 
hotter regions may have altered the initial Be abundance. The 
abundances of Li will help in this. 

Most stars in our sample have been analyzed for Li abun- 
dances before. Our interest in Li in this work is only related to 
the information it might give on the possible effects of convec- 
tive mixing, not on the absolute value of the Li abundance. Thus, 
we decided to adopt Li values from the literature whenever it was 
available. As main reference for lithium abundances (32 stars), 
we adopted the work by Charbonnel & Primas (2005), who an- 
alyzed a large number of previous literature results in detail. 

^ The spectrum is available for download at the ESO website: 
www.eso.org/observing/dfo/quality/UVES/pipeline/solar_spectrum 
.html 
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Fig. 4. Fit to the Li i doublet at 6707.8 A in star HIP 10449. The 
solid line is the observed spectrum and the dashed line the syn- 
thetic one. 



Whenever possible, we adopted their so-called moif consistent 
values including NLTE corrections. As a second preferred refer- 
ence (24 stars), we adopted the results of the survey of Li abun- 
dances by Chen et al. (2001). These Li abundances also include 
corrections for NLTE effects. For stars not included in these pa- 
pers, the Li abundance was obtained from a number of other 
papers: for seven stars, abundances were taken from Boesgaard 
et al. (2005); for four stars, abundances were taken from FOO; 
for three stars, from Romano et al. ( 1999) ; while Favata et al. 
( (19961), G ratton et al. ( i2000l l, Ryan & De hyann is ( [19951 ), Spite et 
al. ('1994), and Takeda & Kawanomoto ( |2005l l each contributed 
with one star. 

Lithium abundances were not available in the literature for 
15 stars of our sample. For these stars, we present new abun- 
dances derived using spectral synthesis. We adopted the atomic 
data for the Li i doublet and neighboring lines listed by Andersen 
et al. ( [T984b . The Li abundance of star HIP 92781 could not 
be calculated because of a damaged profile. The lithium abun- 
dances are listed in table [3] along with a flag indicating the 
sources. An example of a synthetic fit to the Li i doublet is shown 
in Fig. 

Since Li and Be are destroyed at different depths, even stars 
depleted in Li are expected to have preserved their original Be 
abundance (as is likely in the case of the Sun and of cool main 
sequence stars, cfr. e.g. Randich 2007). Whatever the Li refer- 
ence adopted, non-mixed metal-poor dwarfs are expected to have 
a lithium abundance close to the primordial value of the Spite 
plateau (Spite & Spite |1982l l. The exact value of the plateau and 
possible dependencies with [Fe/H] and/or Teff are still controver- 
sial, but for our purpose it is sufficient to assume a typical value 
of A(Li) = 2.20. 

4.4. Aipha elements 

Abundances for the a-elements were taken from the work of 
Venn et al. (I2004I I. These are average abundances of Mg, Ca, and 
Ti determined in the same papers, when available, used as the 
source of the atmospheric parameters we adopt in this work. The 
abundances have not been homogenized since different works 
use different spectral lines to calculate the abundances. We refer 
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the reader to Venn et al. (|2004| l and the original papers for more 
details. We list the abundances in Table |6l We plan to improve 
this analysis by evaluating CNO abundances for all the stars in a 
subsequent work (Smiljanic et al. 2009, in preparation). 

4.5. Uncertainties 
Beryllium 

The Be abundance determination is affected by the uncertainty 
in the atmospheric parameters and uncertainty in the determina- 
tion of the pseudo-continuum. To estimate the effect of the atmo- 
spheric parameters, we change each one by its own error, keep- 
ing the other ones with the original adopted values, and recalcu- 
late the abundances. Thus, we measure what effect the variation 
of one parameter has in the abundances. These effects are listed 
in Table]?] The calculations were done for three stars representa- 
tive of the range of parameters defined by our sample. The total 
uncertainty also includes the effect of the continuum uncertain- 
ties (+ 0.05 dex) and the influence of the unidentified blends 
(+ 0.08 dex) on the Be abundance. We recall that we adopted 
larger uncertainties for Tgff and [Fe/H] when compared to the 
values adopted by FOO and Ed93. Thus, the uncertainties caused 
by these parameters in the abundances already consider some of 
the effects of putting the two samples together 

The uncertainty due to the continuum was determined by 
estimating the sensitivity of the derived Be abundance on the 
continuum level adopted. This uncertainty is mostly related to 
the S/N of the spectrum. After a number of tests with different 
choices for the continuum level, we found this effect to introduce 
a typical uncertainty of + 0.05 dex. 

Another factor that may influence the determination of Be 
abundances are the blends affecting the Be lines. In particular, 
we adopt the value given only by the 3 1 3 1 A line as the abun- 
dance. In this case, the predicted line of Fe i is the most impor- 
tant blend. The influence of this predicted line decreases with 
increasing temperature and/or decreasing metallicity. In a series 
of tests, we noticed that removing the Fe line would increase the 
beryllium abundance in a star with Teff ~ 5200 K and metallicity 
~ -0.90 by 0.15 dex. In a star with Teff ~ 6200 K and metaflic- 

ity 0.70, no effect is seen at all. We adopt a value of + 0.08 

dex, the average between the maximum and minimum effect of 
neglecting blend, to represent this source of uncertainty and ex- 
tend this value to all the sample stars. If our tentative identifica- 
tion of the line is wrong, we are wrongly modeling the effects 
of its variation with [Fe/H] and Teff. This effect, however, is not 
likely to be as strong as the complete neglect of the line. 

Assuming that the effects of the uncertainties from the pa- 
rameters, the continuum, and the blends are independent, we 
may add them in quadrature to estimate the total uncertainty, 
which is also listed in Table |7] 

Lithium 

In the case of lithium, the main uncertainties also come from 
the uncertainties in the atmospheric parameters and from the de- 
termination of the pseudo-continuum. These uncertainties were 
determined in a similar way to the one described above for Be 
and are listed in Table |8] While the main parameter affecting 
the abundance is log g for Be, for Li it is Teff. The uncertain- 
ties we derive are only valid for the Li abundances determined 
in this work, for the values adopted from the literature we refer 
the reader to the original papers. 
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Fig. 5. (a) Diagram showing the difference between the Be abun- 
dances determined in this work and the ones from the literature, 
as a function of the difference in log g. The correlation coeffi- 
cient, p = 0.84, including the outlier, and p = 0.90, excluding it, 
indicates the reality of the linear correlation apparent from the 
figure. The dashed line indicates a linear fit, excluding the out- 
lier, and the solid line is the line x = y. (b) The difference between 
the Be abundances as a function of the difference in [Fe/H]. The 
correlation coefficient, p - 0.24, indicates the lack of correlation 
as apparent from the figure. 

Alpha elements 

As in Venn et al. ( |2004| l. we do not conduct a detailed error anal- 
ysis on the a-elements since they were determined from a vari- 
ety of indicators. We decided to adopt the same representative 
uncertainty as adopted by Venn et al. (,2004j . o-[a/Fe] = ± 0.15 
dex. 

5. Comparison with previous results 

A number of the stars analyzed in this work have beryllium abun- 
dances previously determined in the literature. A detailed com- 
parison of our results with those obtained previously for each 
star is shown in Appendix |A] for the interested reader. As is 
clear from this comparison, there is also a difference in log g 
in most cases where our abundances do not agree with previ- 
ous results. This shows that to establish a consistent and reliable 
gravity scale is important for properly determining the Be abun- 
dance. We have confidence in our results given the good agree- 
ment found with log g derived using the Hipparcos parallaxes. 

In the following, we discuss both the spread in the atmo- 
spheric parameters between the results found in these papers (in 
particular the one of [Fe/H], that might have some influence in 
the spread of the relation between Be and Fe) and the spread of 
the Be abundances. The spread of Be abundances might show to 
what extent the absolute scale of the abundances is reliable. 

5. 1 . Tlie spread of tine parameters 

Although in the comparison shown in Appendix lAl we show re- 
sults from all papers that previously determined Be abundances 
for our sample stars, in what follows we consider only papers 
published after 1997 (including). Early works most likely em- 
ployed data and techniques of lower quality then modem anal- 
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Table 7. The uncertainties of the Be abundance introduced by the uncertainties of the atmospheric parameters. 



Star a Xeff Cr [„„„ g-^ a y^^m] a com. O" bleml O" total 

HIP 53070 + 0.03 ±0.06 +0.01 ±0.00 ±0.05 ±0.08 ±0.12 

HIP 104660 ±0.03 ±0.08 ±0.01 ±0.02 ±0.05 ±0.08 ±0.13 

HD 159307 ±0.04 ±0.08 ±0.01 ±0.03 ±0.05 ±0.08 ±0.13 



Table 8. The uncertainties of the Li abundances, calculated in this work, introduced by the uncertainties of the atmospheric param- 
eters. The total uncertainty also includes the effect of the continuum uncertainties (+ 0.05 dex) on the Li abundance. 



Star g" Teff O" logg 0"f g" [Fe/H] O" com. O" total 

HIP 74067 ±0.10 ±0.00 ±0.00 ±0.00 ±0.05 ±0.11 
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Fig. 6. (a) Diagram showing the difference between the Be abun- 
dances determined in this work and the ones from the literature, 
as a function of the log g values. The correlation coefficient, p 
- 0.14, indicates the lack of correlation also apparent from the 
figure, (b) The difference between the Be abundances as a func- 
tion of the [Fe/H] values. The correlation coefficient, p - -0.26, 
indicates the lack of correlation also apparent from the figure, 
(c) The difference between the Be abundances as a function of 
the log(Be/H) values. The correlation coefficient, p = -0.20, in- 
dicates the lack of correlation also apparent from the figure. 

yses and therefore are prone to greater uncertainties. We also 
do not consider upper limits for this comparison, resulting in 25 
possible comparisons. Multiple comparisons for a given star are 
possible. 

The average difference between our log g values and the ones 
adopted in the literature is 0.00 + 0.28. If we exclude the com- 
parison of our log g of HIP 1 14962, 4.30, with the one adopted 
by Molaro et al. Cl997aJ , 3.44, the average is -0.03 + 0.23. This 
shows that, although no systematic effect seems to be present, a 
wide spread of the log g values adopted in the literature exists. 

Regarding the metallicity, [Fe/H], the average difference be- 
tween our values and the ones adopted in the literature is -1-0.12 
+ 0. 19. If again we exclude the comparison with of HIP 1 14962 
by Molaro et al. ( |1997a| l, the average is H-0. 13 + 0. 17. This com- 
parison seems to indicate a weak systematic effect in the sense 
of our adopted [Fe/H] values to be higher then the ones adopted 
in other analyses. It also indicates the existence of a spread of 
the [Fe/H] values adopted in the literature. This spread is a mea- 



sure of the extent to which the determinations can be considered 
reliable. It will be important to keep this in mind when we dis- 
cuss the relation between log(Be/H) and [Fe/H] in the sections 
to follow. We note that the the uncertainty we assume for the 
[Fe/H] values is similar to the magnitude of the spread found in 
this comparison. 

5.2. The spread on the Be abundances 

The average difference between our Be abundances and the ones 
determined in the literature is -0.06 + 0.19. This indicates that, 
on average, our values are slightly lower than the ones previously 
derived in the literature. As shown in Fig.|5] the difference in the 
Be abundance correlates very well with the difference in log g 
but not with the difference in [Fe/H]. This comparison shows 
the extent to which different analyses of Be can be trusted when 
compared to each other It also shows that most of this spread 
is related to the difference in the adopted log g. Also, no corre- 
lation between the difference in the Be abundance and the own 
abundance, [Fe/H], and log g was found (Fig.|6]l. Therefore, this 
difference does not seem to be due to some systematic effect of 
the analysis. As shown before, our log g values are in excellent 
agreement to the ones derived using Hipparcos parallaxes, giv- 
ing confidence in the abundances we derived. 

6. Be depleted stars 

The light elements Li, Be, and B are very fragile and are de- 
stroyed when in regions inside stars with temperatures of 2.5 x 
10'' K, 3.5 X 10'' K, and 5 x lO*" K, respectively. Main-sequence 
stars with M < 1 .5 M© have a surface convective zone that might 
reach regions hot enough to deplete the surface abundance of 
these elements. The depth of the convective zone is a function of 
effective temperature and metallicity; it is larger when the star is 
cooler and more metal-rich. 

In this section we use Li and Be abundances to evaluate 
which stars are affected by some kind of mixing event. The stars 
that do not display the original Be abundances in their photo- 
spheres will be excluded from further discussion. We also ex- 
clude star HIP 59490, the Be-rich star discussed in detail in 
Smiljanic et al. ( I2008I I. Both Li and Be abundances of this star 
are not a result of the normal evolution of these elements in the 
Galaxy, but a likely result of a peculiar event tentatively identi- 
fied as a hypernova. 

Seven stars in our sample have only Be upper limits. All have 
also low Li abundances. Three of them are subgiants, HIP 17001, 
HIP 71458, and HIP 77946, and have probably diluted both Li 
and Be due to the deepening of the convective envelope, as ex- 
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Fig. 7. Abundances of Li as a function of temperature for the 
sample stars. Detections are shown as filled circles and upper 
limits as upside down open triangles. 



Fig. 8. Abundances of Be relative to Fe as a function of log g for 
the sample stars. A linear fit to the points is also shown. 



pected from stellar evolution. One is HIP 55022 (HD 97916), a 
star first found to be Li depleted by Spite et al. ( 1984). It is one of 
the so-called "ultra-lithium-deficient" stars (Ryan et al. 2001 ), a 
group of a few metal-poor stars that deviate from the Li plateau 
and only have Li upper limits. Ryan et al. ( 1200 11 1 suggest they 
are formed by the same mechanism as forms field blue strag- 
glers. Boesgaard (2007) find these stars are also Be-depleted, 
a result we confirm for HIP 55022. The thick-disk dwarf HIP 
59750 is more metal rich ([Fe/H] - -0.60) than the stars that 
usually define the Li plateau. At its temperature (6200 K), no 
Li (or Be) depletion is expected. Star HIP 36818 has a lower 
temperature (5672 K) and higher metallicity ([Fe/H] = -0.83) 
than plateau stars. Stars of similar temperature and metallicity 
also show Li depletion but normal Be. Both seem to be affected 
by stronger mixing than other similar stars. They might as well 
be metal-rich counterparts of the "ultra-lithium-deficient" stars 
mentioned above. Star HIP 19814 has a rather high Be upper 
limit, consistent with stars of similar temperature and metallic- 
ity where Be was detected. Thus, strictly speaking it cannot be 
considered to be Be-depleted. All these 7 stars will not be con- 
sidered in the discussion presented in the following sections. We 
thus start with a sample of 75 stars where Be was detected. 

Spite & Spite ( |1982| l were the first to show halo stars with 
a range of metallicity to have a constant Li abundance (the so- 
called Spite plateau). Thick-disk, metal-poor dwarfs also show 
the same feature (Molaro et al. |1997bl l. Although departures 
from the plateau due to lithium production by various stellar 
sources seem to start at [Fe/H] ~ -1.00, a sudden large in- 
crease with stars reaching A(Li) ~ 3.0 is observed only for 
[Fe/H] > -0.45, out of the explored range (Travaglio et al. 
1200 11 1. Although some discussion exists on whether the plateau 
is primordial in origin or a result of some depletion mechanism 
(Bonifacio & Mo laro [19971 Piau et al. [2006; Bonifacio et al. 
120071 Korn et al. 120071) . we assume in this work that plateau 
stars did not suffer a previous significant depletion of Li (and 
Be) in their atmospheres. 

In Fig.|2|we show the lithium abundances as a function of 
Teft. As previously found in the literature, stars with Teif > 5800 
K define the Li plateau, while stars with Tefj < 5800 K start to 
show the effects of Li depletion. There are only two stars with 



Teff above 6000 K and A(Li) below the plateau. HIP 85963 is 
a thin-disk star with [Fe/H] = -0.71, Teff = 6227, and A(Li) < 
0.70, and HIP 107975 is a thin-disk star with [Fe/H] = -0.50, 
Teff = 6275, and A(Li) < 1.10. Both seem to be lithium-diifl 
stars or to be evolving from the lithium dip. Since they have 
most likely experienced strong mixing, we do not consider them 
in the further analysis, reducing our sample to 73 stars. Similar 
plots of A(Li) as a function of log g and [Fe/H] do not reveal any 
other star showing strong mixing. 

No trends or deviating stars are seen in plots of [Be/Fe£| as a 
function of Teff and [Fe/H]. We use [Be/Fe] in this case because 
log(Be/H) is expected to vary from star to star, while [Be/Fe] 
is expected to be roughly constant. A small trend with a cor- 
relation coefficient of p = 0.19 is seen in a plot of [Be/Fe] as 
a function of log g (Fig. [8]), suggesting a small decrease in Be 
abundance with decreasing log g. Stars with log g > 4.00 seem to 
form a plateau, while some stars with log g < 4.00 seem to have 
low [Be/Fe]. A similar effect is not seen in the Li abundances, 
which argues against a simple dilution or mixing effect. These 
stars span almost the whole range of metallicity (from -0.50 to 
-1.80), which argues against a chemical evolution effect. The 
trend might be related to the dependence of Be on gravity, sug- 
gesting that log g for the lower gravity stars is possibly not well 
constrained. Nevertheless, we see no strong reason to exclude 
these stars from the discussion. 

Thus, only two stars are excluded from the sample because 
of possible mixing effects. The remaining 73 stars retain their 
original Be abundances and are used in the following sections to 
discuss the galactic evolution of this element and the use of Be 
as a time scale tracing the different star formation history of the 
halo and the thick disk. 

Finally, in Fig. [9] we plot the [Be/Fe] as a function of the 
lithium abundances. No clear trend between Be and Li is seen. 



The so-called lithium-dip is a strong decrease in Li abundances of 
main sequence stars in an interval of ~ 300 K around ~ 6700 K, an 
effect first noticed by Wallerstein et al. (1965 ) in the Hyades and later 
confirmed by Boesgaard & Tripicco ( 19861. A discussion of the physi- 
cal mechanisms possibly associated to the lithium dip can be found in 
Charbonnel & Talon ( 1999 and references therein). 

We adopt the meteoritic solar abundance, A(Be)o = 1.41 (Lodders 
i2003j to calculate the [Be/Fe] ratio. 
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Fig. 9. The [Be/Fe] ratio as a function of the Li abundance. Stars 
with beryllium abundances but lithium upper limits are shown as 
filled circles and stars with both beryllium and lithium detections 
are shown as open circles. A linear fit to the points, without the 
excluded stars, is also shown. No coiTelation between the abun- 
dances of Li and Be is seen. 



The stars with lower Be abundance are not the ones with smaller 
Li abundance. Actually, most of the stars with Li upper-limits 
have Be abundances compatible with stars that show no signs of 
Li depletion. We note that three stars for which we determined 
only a Be upper limit have Li detection, however, they are cool, 
TefF < 5500 K and Li-depleted relative to the plateau. This is 
probably a result of the much higher S/N obtained in the red re- 
gions of the spectra. These stars are nevertheless very interesting 
and should be further investigated. 

7. Galactic evolution of Be 

Once we have cleaned our original sample by identifying the 
stars that have suffered Be dilution or depletion, or for which we 
cannot determine Be abundances, we can investigate the evolu- 
tion of Be in the Galaxy. We are left with 73 stars, 6 from the thin 
disk, 27 from the thick disk, and 39 from the halo. One star has a 
50% probability of being either halo or thick disk. We keep this 
star and the thin disk stars in the sample because according to 
different kinematic criteria they receive a different classification 
(such as in the dissipative and accretion components defined by 
Gratton et al. 2003a - see below). In Fig.[TO]we show a diagram 
of the [Fe/H] vs. log(Be/H) and in Fig. [TT]a diagram of [a/H] 
vs. log(Be/H) for all the 73 stars. In both figures, the points are 
distributed along linear relations, as previously found in the lit- 
erature. We fitted a straight line to the data, taking the errors on 
both axes into account. Following the discussion in Sect. |5] we 
adopted an error of 0.15 dex on [Fe/H] and [a/H]. The adopted 
error on the Be abundance was 0.13 dex. The resulting fits are 
the following: 

log(Be/H) = (-10.38 ±0.08) + (1.24 + 0.07) [Fe/H] (1) 

log(Be/H) = (-10.62 ±0.07) + (1.36 ± 0.08) [q-/H]. (2) 

For the fit with [Fe/H] the goodness-of-fit provides a probability 
of 0. 14, while for the one with [a/H] this is 0.76. Although both 
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Fig. 10. Diagram of [Fe/H] vs. log(Be/H) for all the 73 stars 
showing a linear fit to the points. A correlation coefficient of 
p = 0.89 is found. An example of error bar is shown in the lower 
right corner. 



fits are formally acceptable, the difference in goodness-of-fit be- 
tween the two is remarkable. The probability is very sensitive to 
the eiTor estimates. If we decrease our eiTor estimate on [Fe/H] 
and [a/H] to 0.12 dex, the probability of the two fits drops to 
0.002 and 0.14 respectively. Therefore, while the linear fit with 
[a/H] remains acceptable, the one with [Fe/H] would not, gener- 
ally, be considered acceptable. Visual inspection of Figs. [TOl and 
[TT]confirms the impression derived from the statistical analysis. 
While in Fig.fTTIone can see a clear trend, albeit with an obvious 
scatter, one has the impression in Fig[TO]of being able to see two 
parallel linear relations. 

The most important relation to study the chemical evolution 
of Be is the one with [O/H], since oxygen gives the highest con- 
tribution to the spallation process responsible for forming Be. 
In the absence of oxygen abundances for our complete sample, 
we decided to use the average abundance of a-elements, as com- 
piled by Venn et al. ( |2004| l. These are typically an average of Mg 
and Ca; when available Ti also enters into the definition of a by 
Venn et al. (2004). These elements show behavior with metal- 
licity similar to oxygen. It is therefore reasonable to use [a/H] 
as a proxy for [O/H]. No major inconsistency should be intro- 
duced by this choice. That a single linear relation links Be and 
a elements, even when considering stars extracted from differ- 
ent Galactic populations (halo, thick disk, thin disk) is a conse- 
quence of the uniqueness of the physical process leading to Be 
formation. A non-uniqueness of the relation between Be and Fe, 
if present, could be a consequence of the different evolution of 
Fe in the different Galactic populations. 

In Table |9] we compare the linear relations we derived with 
the ones previously obtained in the literature for metal-poor 
stars. The relation with oxygen or a-elements shows a sizable 
scatter in the fitted slope; this is likely due to the difficulty in 
the measurement of oxygen abundances in metal-poor stars. Our 
own slope is close to the ones obtained by King (2002) using 
[Ca/H] and Boesgaard et al. ( 11999 ) using [O/H]. The zero points 
seem to agree better, the exception being the one determined by 
King ( 2002 ) using [Ca/H] . Our own value is the same as obtained 
by King (2001 ) and Molaro et al. ( |1997al ) both using [O/H]. The 
slope with the a-elements abundance implies the need for a pri- 
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Table 9. A comparison of the linear relations between Be and Fe and O (or any other a indicators adopted) as found in the literature. 
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(1) A mean oxygen abundance from different indicators, such as the molecular OH UV lines and the A 6300 [OI] forbidden line. 
(2) The fits were recalculated in this work based on the original published data. 
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Fig. 11. Diagram of [a/H] vs. log(Be/H) for all the 73 stars show- 
ing a linear fit to the points. A correlation coefficient of p = 0.91 
is found. An example of error bar is shown in the lower right 
corner. 



mary production of Be. All the relations with Fe have a slope 
very close to one. The zero points of the fits have a larger scatter 
The zero point of our fit is closer to the one obtained by King 
(120011 ). who analyzed results available in the literature. 

The possible presence of a real scatter in the Be-Fe relation- 
ship needs to be considered further. The spread is also present 
in earlier data samples but, the limited size of the samples and 
the large errors on Be abundance did not allow a robust assess- 
ment of its existence. In particular, the presence of this scatter 
in samples of smaller size might have motivated the claim of a 
change in the slope, as in Molaro et al. ( |1997al ) and Boesgaard et 
al. ( 119991 ). Molaro et al. ( |1997al l were the first to suggest a pos- 



sible change in the slope of the relation of Be with Fe at [Fe/H] 

1.6, -1.1. No such change is apparent in our data. 

In the hypothesis that Be is a good cosmochronometer, there 
could be a natural explanation for the scatter At a given time 
in the early-Galaxy, the Galactic halo has a higher metallicity 
than the thick disk. Thus, a star formed at this time in the halo 
will have a higher [Fe/H] than a star formed at the same time in 
the thick disk, in spite of a similar Be abundance. We devote the 
next sections to a detailed investigation of the presence of the 
scatter and to its possible causes. This investigation will help to 
constrain the limits of validity of Be as a cosmochronometer. 

7.1. Scatter in the Be-Fe relation 

In Fig. [12] we display our working sample of 73 stars, using 
different symbols to identify several subsets. In addition in the 
plot the two stars, H P 940 28 and HD 132475, identified by 
Boesgaard & Novicki (120061 ) to lie at about 4cr above the mean 
Be-Fe trend are also shown. In this sense these two stars must 
be considered Be-rich. However, we note that HD 94028 is also 
included in our sample. While Boesgaard & Novicki (2006) find 
log(Be/H) = -11.49, we derived log(Be/H) = -11.80 and thus 
do not confirm the 'Be rich status' of this star Although the re- 
sults basically agree within 1 cr (we estimate + 0.13 dex for our 
results and Boesgaard & Novicki + 0.12 for theirs), the differ- 
ence is associated to the different log g values adopted, log g = 
4.20 by us and log g - 4.44 by Boesgaard & Novicki. The log g 
derived using the Hipparcos parallax is log g - 4.27. As we point 
out in Sect. |5] correct determination of log g is essential for de- 
riving consistent Be abundances. In Fig.[T3]we show our fit to the 
spectrum of this star and a synthetic spectrum calculated with the 
larger abundance found by Boesgaard & Novicki (but calculated 
with our adopted atmospheric parameters). Another, even more 
striking exampl e of B e-rich star is HD 106038 (=HIP 59490, 
Smiljanic et al. [2008), which we have removed from our ana- 
lyzed sample. Boesgaard & Novicki (2006) suggests that large 
Be abundance arises from local enrichment phenomena, such as 
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Fig. 12. Our sample of 73 stars in the [Fe/H] - log(Be/H) plane 
split in several sub samples. Circles are thick-disk stars with 
A(Li) > 2.0, squares halo stars with A(Li) > 2.0, x symbols 
are thick-disk stars with A(Li) < 2.0 and + symbols are halo 
stars with A(Li) < 2.0. Asterisks are thin-disk stars or stars with 
ambiguous kinematic classification. The two star symbols are 
the Be-ri ch star s HD 94028 and HD 132475 from Boesgaard & 
Novicki ( [2006] l. HD 94028 (HIP 53070) is also included in our 
sample, where a smaller Be abundance was derived. A short line 
connects the two values. The solid line is a least squares fit to the 
open square symbols and the dashed line is a least squares fit to 
the open circles symbols. A representative error bar is shown at 
the bottom right of the plot. 



may happen in "Super Bubbles" in the vicinity of a Supernova 
(see Parizot l2000. Parizot & Drury '2000' and references therein). 
For the large Be enhancement of HD 106038, Smiljanic et al. 
( 1200 8') found such a scenario insufficient and invoked a hyper- 
nova. Such Be-rich stars seem to be quite rare, and at present it 
seems unlikely that they contribute to a large scatter in Be abun- 
dances. The Galactic orbits of all three stars have been computed 
by Caffau et al. (2005 ) and while HD 132475 and HD 106038 
are halo stars, HD 94028 displays a disk-like kinematics. This 
suggests that the events giving rise to Be-rich stars are indeed 
quite local and do not depend on the stellar population to which 
the star belongs. 

Inspection of Fig.[T2]gives the visual impression that squares 
and circles (halo and thick disk) define two distinct linear re- 
lations. At the same time it is clear that the Li-depleted stars 
would blur this picture, especially those belonging to the halo. 
This is confirmed by fitting a straight line to the halo stars alone. 
The fit to all the halo stars (39 stars) has a probability of 0.04, 
marginally acceptable, but suspicious. If we remove from the 
sample all the Li-depleted stars we are left with a sample of 20 
stars, and the probability of the fit jumps up to 0.19. This fit is 
shown in Fig. [T2]corresponding to 

log(Be/H) = (-10.76 + 0.15) + (0.97 + 0.10) [Fe/H]. (3) 

The sample of thick-disk stars, on the other hand, shows a 
reasonable fit (0.69 probability) without the need to remove the 
Li-depleted stars. The fits obtained by keeping all the thick-disk 
stars and by removing the Li-depleted stars are very similar. For 



Fig. 13. Fit to the region of the beryllium lines in star HIP 53070 
(HD 94028). The crosses represent the observed spectrum, the 
solid line the best fit, and the dashed line a synthetic spec- 
trum calculated with our adopted atmospheric parameters and 
the beryllium abundance determined by Boesgaard & Novicki 
( 12006b . log(Be/H) = -11.49. 

consistency in Fig. [12] we give the one without Li-depleted stars, 
which corresponds to 

log(Be/H) = (-10.30 + 0.25) + (1.21 + 0.20) [Fe/H]. (4) 

The choice of eliminating the halo Li-depleted stars from the 
fit is somewhat arbitrary and we have no particular justification 
for it. Even more so, if we note that many of the halo Li-depleted 
stars have Be abundances that are higher than the mean trend de- 
fined by the other stars. It is nevertheless remarkable that making 
this choice the goodness-of-fit increases. 

Our conclusion is that there is marginal evidence of scatter in 
the Be-Fe relation, above what can be justified by observational 
errors. One possibility is that two distinct Be-Fe relations exist: 
one for the halo and one for the thick disk. Another possibility is 
that there is simply a dispersion in Be abundances at any given 
metal licity. 

7.1 .1 . Stars with similar parameters 

We are analyzing a set of high quality data of stars that have, for 
the vast majority, parallaxes measured by Hipparcos and have 
been quite extensively studied in the literature. Nevertheless, we 
are aware that the uncertainties in an abundance analysis are of- 
ten underestimated. 

We believe a meaningful test of the existence of real scat- 
ter is to identify stars with similar metallicity and similar atmo- 
spheric parameters, hence almost identical spectra, but different 
Be abundances. The difference in strength of the Be lines should 
be obvious and largely independent of our ability to model the 
spectra. 

We identified 4 pairs of stars with similar atmospheric pa- 
rameters and metallicities. Their spectra are compared in Fig. 
[I4l In this figure, all the pairs of spectra show very good agree- 
ment in the behavior of the neighboring metal lines, but clearly 
differ in the two Be lines. This difference is a definitive and con- 
vincing test that the scatter in the abundances is a real feature. 
The stars being compared range from [Fe/H] - -0.80 to -1.40, 
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Fig. 14. Comparison in the Be region between spectra of stars with similar metallicity but different Be abundances, (a) G 5-19 with 
Teff/log g/[Fe/H]/log(Be/H) = 5942/4.24/-1.10/-11.98, shown as a dashed line, vs. HD 199289 with Teff/log g/[Fe/H]/log(Be/H) 
= 5894/4.38/- 1.03/- 11. 36, shown as a solid line, (b) HIP 7459 with Teff/log g/[Fe/H]/log(Be/H) = 5909/4.46/- 1.15/- 11. 90, 
dashed hne, vs. HIP 58145 with Teft/log g/[Fe/H]/log(Be/H) = 5946/4.41/-1.04/-1 1.31, solid line, (c) HIP 58962 with Teff/log 
g/[Fe/H]/log(Be/H) = 583 1/4.36/-0.80/- 11.68, dashed line, vs. G 5-40 with Teff/log g/[Fe/H]/log(Be/H) = 5863/4.24/-0.83/- 11.15, 
sohd line, (d) HIP 22632 with Teff/log g/[Fe/H]/log(Be/H) = 5825/4.30/-1.40/-12.17, dashed hne, vs. HIP 53070 with Teff/log 
g/[Fe/H]/log(Be/H) = 5900/4.20 /-1. 40/- 11. 80, solid line. 



showing that the scatter is not concentrated in a single metal- 7.1.2. Binaries 
licity value, the stars in each pair have very similar gravities, as 
well as similar Li abundances. We cannot therefore invoke differ- 
ences in gravity or in dilution/depletion to explain the observed 
Be difference. 



The situation for a elements is less clear. In the above pairs 
of stars, the [a/H] ratio is not the same, so we could expect stars 
with the same metallicity but different oxygen abundance to have 
different Be abundances. It is more difficult to identify stars in 
our sample with similar atmospheric parameters and also the 
same [a/H], but different Be. Of the stars in Fig.[Tl HIP 22632 
and HIP 53070 are the ones with the the most similar [ff/H], 
-1.04 and -0.96 dex, respectively. In Fig. [15] we compare two 
calcium lines, A 6166.44 A and A 6169.04 A of these stars. As is 
clear from this figure, there is no difference in the calcium abun- 
dances of these stars. The difference in a-element abundances 
of the other stars in Fig. [T4]are on the order of -0. 15, which is 
clodse to the Icr uncertainty of the abundances. The case of HIP 
22632 and HIP 53070 shows that stars with similar abundances 
of a-elements may have different abundances of Be. 



We have finally investigated some other possible factors that 
may influence the scatter, such as the use of data from different 
sources or the possible presence of binaries. We collected infor- 
mation on the multiplicity of our sample stars in the literature, 
which is listed in Table [T] where whenever possible we list the 
number of components of the system (AB for two components, 
ABC for three, and so on). When the number of stars in the sys- 
tem was not clear we only list a flag yes to indicate the star is 
a binary or multiple system. The references consulted employed 
a variety of methods, for example, radial velocities, photometry, 
or speckle interferometry. A flag no does not necessarily mean 
the star is single, but that a reference was found where at least 
one method failed to detect multiplicity. For some stars no refer- 
ence reporting an investigation on multiplicity was found, in this 
case no flag is shown in the table. Again, no difference between 
the subsample of possible binaries and the other stars is found. 
Binaries therefore do not seem to be responsible for the observed 
scatter. 
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Fig. 15. Comparison between two calcium lines in the spectra of 
stars HIP 22632 (dashed line) and HIP 53070 (solid line). The 
Be region of these stars is shown in the last panel of the previous 
Figure. The stars have similar metallicity, similar atmospheric 
parameters, and similar abundances of a-elements but different 
Be abundances. 



7.1.3. The FOO subsample 



As discussed in Sect. |3.2| thei'e is a systematic difference in tem- 
perature and [Fe/H] between the results of Edvardsson et al. 
( I1993I I and Fulbright (2000), although within the uncertainties,. 
This difference in metallicity may introduce some scatter into 
the relations shown in Figs. [10] and [TT] In the final sample of 73 
stars, 49 have parameters calculated by FOO, of which 26 are halo 
stars and 19 are thick-disk stars. We repeated the analysis with 
this subsample of stars. The linear fits using only the Fullbright 
stars and rms are statistically identical to the ones obtained with 
the whole clean sample (also listed in Table |9]l. The equivalent 
of Figs.[TO]and[TT]for the FOO subsample are given in Appendix 

IB 



7.2. log(Be/H) vs. [Fe/Be] 

The evolution of a given element in the Galaxy is usually ana- 
lyzed in a plot of [Element/Fe] vs. [Fe/H]. If Be is a better time 
scale than Fe, then changing Fe for Be in this kind of discussion 
would be then the natural follow-up. In Fig. [16] we plot [Fe/Be] 
as a function of log(Be/H). 

A decrease in [Fe/Be] with increasing log(Be/H) is seen. It 
is interesting to see that for log(Be/H) > -11.5 the majority of 
the stars have a sub-solar Fe/Be ratio. In principle similar plots 
could be constructed also for elements other that iron. 



8. Be as a chronometer 

The hypothesis that Be can be used as an age indicator relies 
on assumptions about the uniformity of the GCR in the early- 
Galaxy. This uniformity is necessary to assure the Be clock runs 
at a uniform pace within the relevant region of the Galaxy. In 
general, models of Galactic chemical evolution do not include 
a detailed treatment of GCR propagation and confinement. The 
model we use to compare with our observations (see Valle et 
al. [2002 and references therein) follows a multizone treatment 
of the Galaxy, where each component is characterized by dif- 




Fig. 16. Diagram of log(Be/H) vs. [Fe/Be]. Thick disk stars are 
shown as filled circles, halo stars as starred symbols, and the 
subgroup of 'Fe-rich' halo stars discussed Sect. [8] is shown as 
open circles. 



ferent star formation rates with a radial dependence. The GCR, 
however, has neither radial nor vertical dependence. The inter- 
pretation of the results in light of these models may, therefore, 
be somewhat limited. 

Pasquini et al. (2005) began using the Be abundances as a 
time scale, and show that halo and thick-disk stars separate in a 
diagram of log(Be/H) vs. [a/Fe]. Their sample, however, is com- 
posed of only twenty stars analyzed by Boesgaard et al. (I1999I I. 
In a diagram [O/Fe] vs. log(Be/H), the beryllium abundance in 
the abscissa can be considered as 'increasing time' while the 
[O/Fe] ratio in the ordinate can be considered as the 'star for- 
mation rate' (SFR). We use [a/Fe] here rather than [O/Fe]. In 
Fig. [it] we show the [a/Fe] vs. log(Be/H) for our sample. The 
figure is divided into two panels; in the upper panel we plot the 
thick disk stars and in the lower panel we plot the halo stars. 
In the same figure the two models from Pasquini et al. (120051 1. 
derived from Valle et al. ( 120021 ), are also shown. In Fig. [TSlwe 
show the classical [a/Fe] vs. [Fe/H] diagram, also divided with 
the thick-disk stars in the upper panel and the halo stars in the 
lower panel. 

From Fig.[l7]it appears that the model is fairly successful in 
reproducing the general trend of the thick-disk stars. However 
the halo stars appear to split along a double sequence. The 
most striking is the group of stars with high [a/Fe] and high 
log(Be/H), which the model totally fails to reproduce. From the 
point of view of a elements and beryllium abundance, this group 
of stars is essentially indistinguishable from the thick-disk stars. 

The first conclusion that may be drawn from this plot is that 
the thick disk appears to be a fairly homogeneous population, 
where the "Be chronometer" may be used with some confidence. 
The halo instead appears to be more complex. In the follow- 
ing sections we discuss the use of Be as a chronometer in the 
two populations separately. Galactic orbits for the majority of 
the sample stars were computed by Gratton et al. (,2003 a). For 
the few remaining stars, new orbits were calculated following 
the same method. We refer to Gratton et al. ( |2003al ) for the de- 
tails. The orbital parameters are listed in Table[6] In Figs. [T9l and 
l20l we show the correlation of chemical and dynamical proper- 
ties. Alternative to the halo - thick disk distinction, we could 
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Fig. 17. Diagram of [a/Fe] vs. log(Be/H). The thick-disk stars 
are shown in the upper panel, while the halo stars are shown in 
the lower panel. The thick-disk stars diagram is characterized by 
a large scatter, while the halo stars diagram clearly divides into 
two sequences. The curves are the predictions of the models by 
Valle et al. (|2002] l. 



have used the accretion - dissipative component distinction in- 
troduced by Gratton et al. ( |2003al|2003b] l. We refer to the orig- 
inal papers for a detailed presentation of the kinematic criteria 
used for this division. We did this exercise, and the main con- 
clusions are essentially the same. The corresponding plots are 
provided in AppendixICl 

8.1. The Be chronometer in the thick disk 

Thick-disk stars show smaller scatter in all diagrams, and it 
seems possible to understand their formation in a relatively sim- 
ple way. The most interesting feature is, however, the anticor- 
relation of [a/Fe] with Rmin in our sample (Fig. |20] |. At first 
glance, this anticorrelation seems to be driven by the two stars 
with largest Rmin- However, even when restricting the sample to 
stars with R^in < 6 kpc, the probability of correlation is 99.8%. 
A similar anticorrelation is also present in the data of Gratton et 
al. (I2003bl see their Fig. 6). If Rmi„ is representative of the ra- 
dius of formation of a star, neglecting the phenomena of orbital 
diffusion, then the simplest interpretation is that the thick disk 
experienced a higher SFR in the inner regions than in the outer 
regions. This does not contrast with the notion that the thick is 
very old and massive (Fuhrmann [19981 |2000l |2004l |2008] l. In 
fact the star formation is supposed to have continued for about 2 
Gyr (Fuhrmann 120041) . which is ample time for the external re- 
gions experiencing a low SFR to develop a lower a/Fe than the 
inner regions where the SFR was higher Note also that the most 
Be-poor, i.e. oldest, stars are found at low Rmin- The emerging 
picture is that of a dissipative thick disk formed inside-out. 

Such an anticorrelation is not present in the log(Be/H) vs. 
Rmin plot, as shown in the upper panel of Figure |20| It is also 
not present in the [Be/Fe] vs. Rmin plot. This lack of anticorrela- 
tion is an important result. The major objection to the use of Be 
as a time indicator is the possibility that the abundance of this 
element is dominated by local events rather than by the relativis- 
tic component. The log(Be/H) vs. R^in diagram does not show 
evidence of a general decrease in Be with the distance from the 
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Fig. 18. Diagram of [a/Fe] vs. [Fe/H]. The thick-disk stars are 
shown in the upper panel, while the halo stars are shown in the 
lower panel. The curves are the predictions of the models by 
Valle et al. (l2002li 
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Fig. 19. Diagram of log(Be/H) vs. the V, the component of the 
space velocity of the star in the direction of the disk rotation. 
The thick-disk stars are shown as filled circles, halo stars are 
shown as starred symbols, and the subgroup of halo stars as open 
circles. A typical error of + 12 Km s ' in V was adopted (see 
Gratton et al. i2003bi) . 



galactic center, as could be expected if the local flux of GCR was 
the dominant factor. 

The range in Be covered by the thick-disk stars of our sam- 
ple is relatively small. Moreover, most of the Be-poor stars are 
present at a relatively small Rmin- This agrees with the notion of 
a rather fast formation for the thick disk (Fuhrmann |2004i [20081 1 . 
As a reference time scale, the reader may consider that the thick- 
disk models shown in the upper panel of Fig. [iTj have a star for- 
mation period of 1 Gyr. 
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Fig. 20. Diagram of log(Be/H) vs. R,,,,-,,, the perigalactic distance 
of the stellar orbit. The thick disk stars are shown as filled circles, 
halo stars are shown as staiTed symbols, and the subgroup of halo 
stars as open circles. A typical error of + 0.40 Kpc in Rmin was 
adopted (see Gratton et al. l2003bl l. 

8.2. The Be chronometer in the halo 

For a long time the debate on the halo formation was centered 
on the opposite views of "monolithic collapse" (Eggen et al. 
|1962| l and "fragment assembly" (Searle & Zinn 1978 ). It is now 
understood that these represent rather idealized extreme cases 
that do not happen in reality. Major insights into halo assembly 
have come from simulations of galaxy formation in a cosmolog- 
ical context (see e.g. Font et al. 2006 and references therein). In 
these simulations in the ACDM scenario, the galaxy formation 
proceeds hierarchically, nevertheless they predict the existence 
of metallicity gradients and other features typical of the "mono- 
lithic collapse". In light of these simulations, we should expect to 
find, in a halo sample, stars that have been accreted from many 
dwarf satellite galaxies, as well as stars formed on short time 
scales in the halo (i.e. a "monolithic collapse"-like component). 

Even before considering our data one can expect that the 
"Be chronometer" was running at a different pace in the differ- 
ent components out of which the halo was assembled. This in 
fact appears quite obvious from Fig. [IT] The stars with high Be 
abundances and high [a/Fe] cannot be reproduced by our model. 
It is tempting to interpret these stars as "accreted" from satellites 
in which the Be production was more efficient than in the Milky 
Way. However, these stars do not display any peculiar kinematic 
property. 

The good news is that the halo stars do not fill uniformly the 
log(Be/H) - [a/Fe] plane, but position around only two, rather 
well-defined sequences. This means that there is some prospect 
for also using the "Be chronometer" for the halo stars. Hopefully 
future models will be able to explain the two sequences and cal- 
ibrate the Be abundance versus time. By use of a diagram like 
Fig.[T7] one should be able to decide which sequence any given 
star belongs to. 

The reality of the separation seen in Fig.[T7]is dependent on 
the eiTor of the abundances. A statistical test shows that the halo 
stars in the top right of panel b in this figure have probabilities 
between 40% to 9% of following the model curve. This suggests 
that the separation seen in the figure is real. These probabilities 
are likely to improve, and more stringent constraints might be 



derived, when oxygen, which is a pure SN II product, is deter- 
mined in a homogeneous way for all the sample stars (Smiljanic 
et al., in preparation). 

In Figs. [T9l and |20l we plot both [c/Fe] and log(Be/H) as a 
function of V, the spatial velocity in the direction of the rota- 
tion of Galaxy, and Rmh,, the perigalactic distance, respectively. 
At the end of the '90s, several halo stars with a low [a/Fe] were 
identified (Carney et al. 1997 Nissen & Schuster 1997) and were 
considered as "accreted" objects. Inspection of Fig. [17] shows 
that such stars are however predicted by our evolutionary model. 
To easily identify the "low a" stars (which we define as stars 
with [a/Fe] < 0.25 and log(Be/H) < - 11.4) in the figures, we 
plot them as open circles. These stars have very similar kine- 
matic properties. They tend to have V close to zero and Rmin < 
1 kpc, although the statistics are too few to allow a robust claim. 
They seem to be a group of non-rotating stars going very close 
to the Galactic center, a behavior that might be expected to be 
shown by accreted stars that sink to the Galactic center by dy- 
namical friction. They also show a remarkably narrow range of 
Be abundances, suggesting a narrow range in age for this com- 
ponent. These characteristics might hint at a possible common 
origin for these stars. 

A robust interpretation of these components in light of our 
evolutionary model does not seem possible. Unlike what is as- 
sumed by the model, our results suggest that the halo is not a 
single uniform population where a clear age-metallicity relation 
can be defined. The splitting of the halo into the two components 
identified in this work might be related to the accretion of exter- 
nal components or simply to variations in the star formation in 
different and initially independent regions of the early halo. 

9. Summary 

We have presented the largest sample, to date, of Be abundances 
in stars. We confirm the existence of a linear increase in Be with 
[Fe/H] and [a/H], as has been found by all the previous investi- 
gations. However, thanks to the large dimension of our sample 
and the homogeneous analysis, we have a marginal detection of 
an excess scatter in the log(Be/H) - [Fe/H] relation, above what 
may be expected by observational eiTors. A similar scatter is not 
obvious in the log(Be/H) - [a/H] relation, although we identify 
a pair of stars with similar parameters and similar [a/H] but with 
a 0.4 dex difference in Be abundance. Our interpretation of this 
scatter is that it arises from the different evolutionary time scales 
of different Galactic components (halo, thick disk). Beryllium is 
a good indicator of time, anywhere, but metal enrichment pro- 
ceeds at a different pace in different environments. 

We tested the use of Be as a chronometer by comparing our 
observations to model predictions. The result appears to be fairly 
satisfactory for the thick disk, which we confirm to be a homoge- 
neous population. There appears to be no trend in Be abundance 
with perigalactic distance. We take this as evidence that the Be 
production is dominated by the relativistic component of the cos- 
mic rays and largely insensitive to local events. Our models of 
Be evolution, however, fail to reproduce the observations for a 
subset of halo stars. More sophisticated models that explicitly 
consider the merger history of the halo should be envisaged. In 
a plot of log(Be/H) versus [ff/Fe] (Fig.[T7]), the halo stars do not 
fill the plane randomly, but align in two rather clearly defined se- 
quences. Only one of these two sequences can be explained by 
our models. However, most stars of this same component show 
very similar kinematics and a narrow range in Be (age). They 
seem to be a group of non-rotating stars going very close to the 
Galactic center so might have a common origin. 
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For the halo stars we do not identify any significant corre- 
lation between chemical and kinematic properties. For the thick 
disk, instead, we find a significant anticorrelation of [a/Fe] with 
perigalactic distance. This anticorrelation might be interpreted 
as evidence that the SFR was lower in the outer regions of the 
thick disk, pointing towards an inside-out formation. 

Even though not all our observations can be readily inter- 
preted in terms of existing models, we believe that our study has 
highlighted the usefulness of Be abundance in studying the dif- 
ferent Galactic populations. 
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4.31 


A 1 A 

-0.19 


/I 1 O 

4.12 


1 A O C 

+0.25 


2 

no 


1 /I ^ 




£>JJ-zl 34zU 


1 r\ 1 ^ 
lU. 1 J 




1 A A 
i .44 


J.O/ 


A 1 

— U. 10 


3.00 


A /I T 

U.43 




JoVoz 


1U5UU4 




lU.zl 


2.00 


A AC\ 

4.4y 




AIT 

—0.1 / 






AdC 


CO/1 OA 


lUOUJO 


Cj Iz-zl 


10. lo 


O 1 A 

9. lo 


1 cr\ 
1.5U 


A OO 

4.99 


A 1 

—O.lo 


4.83 


A A'2 

—0.03 


2 

no 


cm c A 


1 A/; CI/; 
106516 


HK 4o5 / 


1 1 
0.1 1 


A A '~1 A 

44.34 


1 A 1 
1.01 


4.34 


A 1 1 

-0.1 1 


/I OT 

4.23 


1 A O 1 

+0.21 


A "D 1 

AB 


60632 


108177 


G 13-35 


9.67 






A on 

4.87 


AIT 

-0.17 


/I T A 

4.70 


1 A AO 

H-0.02 


2 

no 


62882 


111980 




8.38 


12.48 


1.38 


3.86 








ab' 


/C? CCA 

63559 


1 1 OAOO 

1 13083 




O AC 

8.05 


1 O C 1 

18.51 


1 1 O 

1.12 


A T A 

4.39 


A 1 O 

-0.18 


/I O 1 

4.21 


1 A OO 

+0.22 


yes"* 


U J" i o 


1 1 ^679 




9 70 




1 ^7 

i .J A. 


^ 87 


W. i u 






4 

no 


64426 


1 14762 




7.31 


24.65 


1.44 


4.27 


-0.16 


4.11 


H-0.26 


12 

yes 


66665 




BD-H13 2698 


9.37 


7.44 


1.70 


3.73 


-0.20 


3.53 


+0.49 


2 

no 


67655 


120559 




7.97 


40.02 


1.00 


5.98 


-0.26 


5.72 


-0.39 




67863 


121004 




9.04 


16.73 


1.35 


5.16 


-0.16 


5.00 


-0.10 


4 

no 


70681 


126681 




9.33 


19.16 


1.44 


5.74 


-0.22 


5.52 


-0.31 


no4 


71458 


128279 




7.97 


5.96 


1.32 










4 

no 



(1) Dommangel^& Nys ( 1994). (2) Latham et al. (2002). (3) Fouls ( I1987D , ( 4) Nord strom et al. (|20041l, (5) Th is work," 
(6) Lu et al. i\9WT\, (7) Lindgren & Ardeberg ( 1996), (8) Gliese & Jahreiss ( 1988), (9) Goldberg et al. ("SU®, 



(10) Carney et al. ( |200TT), (1 1 ) Stryker et al. a985i ). (12) Patience et al. (2Q02), (13) McAlister et al. (|1987i i. 
(14) MacConnell et al. ( fT997] l 



Table 1. continued. 
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HIP 


HD 


Other 


V 


TT 




My 


EC 


Mboi 


log (L JLq) 


binary? 






name 


mag. 


mas 


mas 


mag. 


mag. 


mag. 






72461 




BD+26 2606 


9.74 


10.28 


1.42 


4.80 


-0.21 


4.59 


-1-0.06 


11 

yes 


74067 


134088 




8.00 


28.29 


1.04 


5.26 


-0.20 


5.06 


-0.12 


4 

no 


74079 


134169 




7.68 


16.80 


1.11 


3.81 


-0.16 


3.65 


H-0.44 


4 

yes 


1 /y46 






o 1 
8.61 


O.JO 


1 oo 
1.23 


o /;a 
2.69 


All 

-0.11 


o c o 
2.58 


1 A OT 
+0.8 / 




oUUUi 




/-I 1 AO /II 


11 CI 


y.iz 


A 1 

3.U1 


A 1 

0.31 


A O 1 

— U.zl 


tZ 1 A 

0. lU 


A C/1 

— U.54 


2 

yes 


oUoj / 


i4ooi0 


1 T O 1 


n oo 


0/1 1.A 

z4o4 


A OA 


/I O 1 

4. /I 


AIT 
— U.l / 


A A/1 

4.U4 


1 A OO 

-I-U.2S 


2 

no 


oil nf\ 
Oil /O 


1 /I n /1 1 /I 
149414 




9.63 


OA n 1 
z(J. / 1 


1 C A 

1.50 


o 1 
6.21 


A o/; 
-0.26 


C AC 

5.95 


A /I O 

-0.48 


ArSCUb 


sjy63 


1 cm AT 




n A ^ 


1 O /I A 


A AA 

0.99 


O AC 

3.05 


A 1 O 

-0.13 


O AO 


1 A TO 

-t-0. /3 


4 

yes 


a IWi 




r>U+20 3603 


A -7 1 
9. /I 


0.4/ 


T O C 

/.85 


o n^z 
3. /6 


A 1 A 

-0.19 


O CO 

3.58 


1 A /I T 

-1-0.4/ 


4 

yes 


O O A 1 A 

soOiO 


1 /CO O 1 A 

loiolU 




A iZI 

9.63 


1 1 oo 
1 1.00 


O O 1 

2.21 


C AA 

5.00 


A O O 

-0.28 


A TO 

4. /2 


1 A A 1 

-1-0.01 


Ars 


OT7Q 1 


1 1 TO 
1 /jl /V 




O AO 


1 1 CC 
1 l.JJ 


1 01 

l.al 


/I /I c 

4.45 


A 1 

— U.lo 


A OO 

4./y 


1 A 1 

H-U. lo 


2 

no 


A /I /I /I A 

y444y 


1 /Vozo 




A 1 -7 
9.1 / 




1.36 


c c 

3.55 


A O 1 

—0.21 


O O /I 

3.34 


1 A c/; 

H-0.56 


2 

no 


AOA'^A 

yoo/0 


1 O O C 1 A 




O O "3 

8.83 


'O C TO 


1 IT 
1.1 / 


c o c 
5.85 


A OT 

—0.2 / 


ceo 
5.58 


A O 

-0.33 


4 

yes 


AO 


1 OACCO 




1 .1 1 


14. /o 


1 1 A 
1.10 


CT 

3.5 / 








A T34 

Ars 


1 AA</^C 


1 moA 1 




/^c 


zz.oo 


1 0/1 

l.z4 


C /I c 

5.45 


A 1 O 

— U. 19 


c o/^ 
5.20 


A OA 

— u.zu 


2 

no 


iuu/yz 




1/1 1 c 


8.3o 


1 / .y4 


1 0/1 

l.z4 


A /^O 

4.03 


AIT 

—U.l / 


A A(Z 

4.4o 


1 A 1 O 


2 

no 


1 A 1 o /I 
101346 


1 AC/COO 




O CO 

8.j3 


o.o3 


1 1 
1.16 


O O 1 

3.21 


A 1 O 

-0.13 


O AO 

3.08 


1 A izn 
+0.6/ 


4 

no 


1 AO /I AO 

10J4ya 


1 AAOOA 




O OA 

8.29 


1 O A/1 

io.y4 


1 AO 

1.03 


4.68 


AIT 

-0.1 / 


/I C 1 

4.51 


1 A 1 A 

+0.10 


4 

no 


1 A/1 

104660 


OA1 OOA 

zUlooV 




O A/1 

8.04 


1 T AC 

1 /.yj 


1 A A 

1.44 


/I O 1 

4.31 


A OA 

-0.20 


/111 
4.1 1 


1 A O/C 

+0.26 


A "D 1 

AB 




zUioUo 


TJD 0101 


/I oo 
4.// 


1 AO CA 


A CO 

U.59 


4.4U 


A 1 

— U. 13 


A OT 
4.Z / 


1 A 1 O 

+U. 19 


4 

no 


1 AC O O O 




O C OA 

(j 25-2y 


A AO 

9.03 


1 O AO 

13. Uz 


111 
1.11 


A f\n 
4.0/ 


A 1 

-0.16 


O A 1 

3.yi 


1 A O O 

+0.33 


2 

no 


1 A/C A 

106447 




(j 26-12 


1 O 1 c 

12.15 


1 CT 

-1.57 


A C A 

4.54 




A O 1 

-0.21 






2 

no 


107y75 


207978 


15 Peg 


5.54 


36.15 


0.69 


3.32 


-0.11 


3.21 


+0.62 


13 

no 


108490 


AO A A^ 

208906 




6.95 


34.12 


0.70 


4.62 


-0.15 


4.47 


+0.1 1 


AB' 


109067 




G 18-28 


9.55 


21.52 


1.59 


6.21 


-0.26 


5.96 


-0.48 


AB^ 


109558 




G 126-62 


9.47 


8.43 


1.42 


4.10 


-0.18 


3.92 


+0.33 


AB^ 


109646 


1 AT 

210752 




7.40 


26.57 


0.85 


AAA 


-0.14 


4.30 


+0.18 


4 

no 


1 12229 


215257 


G 27-44 


7.40 


23.66 


0.97 


AAA 


-0.14 


AAA 


+0.25 


2 

no 


1 1 /I 1 

i 142 / i 


O 1 O C AO 

218502 




O C A 

8.50 


1/1 o o 

14.33 


1 OA 

1.20 


A O O 

4.28 


All 

-0.11 


A T n 

AA 1 


1 A o/; 
+0.26 




11/1 o/^o 


O 1 O/^ 1 T 


U / /3-i 


1 T 
0.1/ 


1 O A/1 

1/.U4 


O /I 1 

Z.41 


CT 

3.5 / 


A OA 

— U.zU 


OT 

3.3 / 


1 A CC 

+U.55 




1 15167 




G 2y-23 


10.21 


3.26 


2.20 


2.78 


-0.20 


2.58 


+0.87 


AB' 


117041 


222766 




10.15 


8.46 


1.76 


4.79 


-0.26 


4.53 


+0.09 


no^ 






G05-19 


11.12 








-0.16 






2 

no 






G05-40 


10.79 








-0.17 






2 

no 






G66-51 


10.63 








-0.26 






2 

no 






G166-37 


12.66 


5.20"' 


0.70 


6.24 


-0.27 


5.97 


-0.49 


2 

no 






G 170-21 


12.51 








-0.21 






2 

no 



(1) Dommanget & Nvs ( 1994). (2) Latham et al. (2002), (3) Fouts ( 1987), (4) Nordstrom et al. ( I2004I I, (5) Th is work^' 
(6) Lu et al. (il987 i), (7) L indgren & Ardeberg ( 1996), (8) Gliese & Jahreiss ( 1988), (9) Goldberg et al. ('SUOl), 
(10) Carney et al. dSUOTI), (1 1 ) Stryker et al. a985j , (12) Patience et al. (i2002j , (13) McAlister et al. a987i ), 
(14) MacConnell et al. ( fT997] ) 
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Table 2. Log book of the observations. 



fflp 


Date of obs. 


TTT 1 

Type' 


Exp. Time 


R 




S/N^ 








(s) 




(A) 


(final) 


171 


13.oct.2005 


New UVES 


2x150 


35 000 


3460 


100 




13.oct.2005 


New UVES 


2x30 




5800 






13.oct.2005 


New UVES 


2x30 




8600 




3026 


08.oct.2001 


Archive 


4200 


35 000 


3460 


130 


7459 


22.sep.2005 


New UVES 


2x1800 


35 000 


3460 


60 




22.sep.2005 


New UVES 


2x800 




5800 






22.sep.2005 


New UVES 


2x800 




8600 




10140 


14.oct.2005 


New UVES 


2x900 


35 000 


3460 


70 




14.oct.2005 


New UVES 


2x350 




5800 






14.oct.2005 


New UVES 


2x350 




8600 




10449 


14.oct.2005 


New UVES 


2x900 


35 000 


3460 


90 




14.oct.2005 


New UVES 


2x350 




5800 






14.oct.2005 


New UVES 


2x350 




8600 




11952 


13.oct.2005 


New UVES 


1500 


35 000 


3460 


75 




13.oct.2005 


New UVES 


2x650 




8600 




13366 


13.oct.2005 


New UVES 


450 and 900 


35 000 


3460 


90 




13.oct.2005 


New UVES 


2x150 




5800 






13.oct.2005 


New UVES 


2x700 




8600 




14086 


23.sep.2005 


New UVES 


2x150 


35 000 


3460 


70 




23.sep.2005 


New UVES 


2x30 




5800 






23.sep.2005 


New UVES 


2x30 




8600 




14594 


26.nov.2001 


Archive 


1800 


40000 


3460 


150 


17001 


13.dec.2001 


Archive 


4x1175 


40 000 


3460 


110 


17147 


26.nov.2001 


Archive 


480 


40000 


3460 


140 


18802 


28.nov.2001 


Archive 


1500 


35 000 


3460 


100 


19007 


16.oct.2005 


New UVES 


2x1200 


35 000 


3460 


50 




16.oct.2005 


New UVES 


2x500 




5800 






16.oct.2005 


New UVES 


2x500 




8600 




19814 


16.oct.2005 


New UVES 


2x2400 


35 000 


3460 


55 




16.oct.2005 


New UVES 


2x1100 




5800 






16.oct.2005 


New UVES 


2x1100 




8600 




21609 


23.sep.2005 


New UVES 


2x1800 


35 000 


3460 


60 




23.sep.2005 


New UVES 


2x800 




5800 






23.sep.2005 


New UVES 


2x800 




8600 




22632 


21.sep.2005 


New UVES 


2x900 


35 000 


3460 


60 




21.sep.2005 


New UVES 


2x350 




5800 






21.sep.2005 


New UVES 


2x350 




8600 




24030 


08.oct.2005 


New UVES 


2x1500 


35 000 


3460 


75 




08.oct.2005 


New UVES 


2x650 




5800 






08.oct.2005 


New UVES 


2x650 




8600 




24316 


22.sep.2005 


New UVES 


2x1200 


35 000 


3460 


70 




22.sep.2005 


New UVES 


2x500 




5800 






22.sep.2005 


New UVES 


2x500 




8600 




31188 


26.nov.2005 


New UVES 


2x700 


35 000 


3460 


80 




26.nov.2005 


New UVES 


2x300 




5800 






26.nov.2005 


New UVES 


2x300 




8600 




31639 


24.sep.2005 


New UVES 


2x1200 


35 000 


3460 


25 




24.sep.2005 


New UVES 


2x500 




5800 






24.sep.2005 


New UVES 


2x500 




8600 




33221 


01.dec.2001 


Archive 


4200 


35 000 


3460 


145 




26.dec.2001 


Archive 


4200 




3460 






27.dec.2001 


Archive 


4200 




3460 




33582 


10.nov.2005 


New UVES 


2x900 


35 000 


3460 


65 




10.nov.2005 


New UVES 


2x350 




5800 






10.nov.2005 


New UVES 


2x350 




8600 




34285 


10.apr.2000 


Archive 


4500 


40000 


3460 


100 




ll.apr.2000 


Archive 


3600 


50000 


3460 




36491 


24.dec.2005 


New UVES 


2x500 


35 000 


3460 


40 



(1) Whether the spectra are from new observations, the ESO archive data, or the UVES POP Ubrary. 

(2) The S/N per pixel of the final combined spectrum in the Be region. 
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Table 2. continued. 



T 1 1 ri 

Hlr 


Date 01 obs. 


Type 


Exp. lime 


K 












(s) 




(A) 


(final) 




24.dec.2005 


New UVES 


2x250 




5800 






24.dec.2005 


New UVES 


2x250 




8600 




36640 


08.oct.2005 


New UVES 


2x150 


35 000 


3460 


130 




08.oct.2005 


New UVES 


2x30 




5800 






08.oct.2005 


New UVES 


2x30 




8600 




36818 


25.nov.2001 


Archive 


5400 


35 000 


3460 


75 


36849 


25.nov.2001 


Archive 


2x2100 


35 000 


3460 


130 


37853 


08.apr.2001 


UVES POP 


2x60 


80 000 


3460 


80 


38625 


20.dec.2005 


New UVES 


4x360 


35 000 


3460 


60 




20.dec.2005 


New UVES 


6x120 




5800 






20.dec.2005 


New UVES 


2x120 




8600 




42592 


20.dec.2005 


New UVES 


2x900 


35 000 


3460 


90 




20.dec.2005 


New UVES 


2x350 




5800 






20.dec.2005 


New UVES 


2x360 




8600 




44075 


14.apr.2001 


UVES POP 


200 


80 000 


3460 


60 


44124 


01.feb.2005 


Archive 


3x1066 


40 000 


3460 


105 




01.feb.2005 


Archive 


3x800 




5800 




45554 


24.dec.2005 


New UVES 


2x3000 


35 000 


3460 


45 




24.dec.2005 


New UVES 


2x1400 




5800 






24.dec.2005 


New UVES 


2x1400 




8600 




48152 


29.nov.2002 


UVES POP 


3x400 


80 000 


3460 


115 


50139 


05.jaii.2006 


New UVES 


2x400 


35 000 


3460 


90 




05.jan.2006 


New UVES 


2x130 




5800 






05.jan.2006 


New UVES 


2x130 




8600 




52771 


20.mar.2005 


Archive 


3x1066 


40 000 


3460 


40 


53070 


05.jan.2006 


New UVES 


600 


35 000 


3460 


60 




05.jan.2006 


New UVES 


2x200 




8600 




55022 


10.apr.2000 


Archive 


3x1800 


40 000 


3460 


90 


57265 


14.jan.2006 


New UVES 


2x2100 


35 000 


3460 


50 




14.jan.2006 


New UVES 


2x900 




5800 






14.jan.2006 


New UVES 


2x900 




8600 




58145 


14.jan.2006 


New UVES 


2x1800 


35 000 


3460 


85 




14.jan.2006 


New UVES 


2x800 




5800 






14.jan.2006 


New UVES 


2x800 




8600 




58962 


12.jan.2002 


Archive 


2x4800 


35 000 


3460 


70 


59490 


12.apr.2000 


Archive 


2x3600 


40 000 


3460 


70 


59750 


30.jan.2006 


New UVES 


2x150 


35 000 


3460 


100 




30.jan.2006 


New UVES 


2x30 




5800 






30.jan.2006 


New UVES 


2x30 




8600 




60632 


28.mar.2004 


Archive 


3x1066 


40000 


3460 


70 


62882 


10.apr.2000 


Archive 


2x1800 


40 000 


3460 


155 


63559 


19.jan.2006 


New UVES 


2x300 


35 000 


3460 


75 




19.jan.2006 


New UVES 


1x200 




5800 






19.jan.2006 


New UVES 


2x200 




8600 




63918 


10.apr.2000 


Archive 


2x3600 


40 000 


3460 


100 


64426 


08.jun.2006 


New UVES 


3x450 


35 000 


3460 


145 




08.jun.2006 


New UVES 


4x150 




5800 






08.jun.2006 


New UVES 


2x150 




8600 




66665 


15.jun.2006 


New UVES 


2x1000 


35 000 


3460 


50 




15.jun.2006 


New UVES 


2x400 




5800 






15.jun.2006 


New UVES 


2x400 




8600 




/ODD 


zi.reb.zUUz 


Archive 


zxyoo 


J J 000 


j400 


1 r\f\ 
lUU 


67863 


ll.apr.2001 


Archive 


3600 


40000 


3460 


145 


70681 


09.apr.2000 


Archive 


2x2700 


50000 


3460 


100 


71458 


04.feb.2002 


Archive 


5x250 


40000 


3460 


170 


72461 


25.jun.2006 


New UVES 


2x1800 


40000 


3460 


60 




25.jun.2006 


New UVES 


2x800 




5800 






25.jun.2006 


New UVES 


2x800 




8600 





(1) Whether the spectra are from new observations, the ESO archive data, or the UVES POP Ubrary. 



(2) The S/N per pixel of the final combined spectrum in the Be region. 
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Table 2. continued. 



XJTT> 

rllr 


Date of obs. 


Type 


Exp. Time 


K 


1 


o/JN 








(s) 




(A) 


(final) 


74067 


12.jun.2006 


New UVES 


2x500 


35 000 


3460 


100 




12.jun.2006 


New UVES 


2x200 




5800 






12.jun.2006 


New UVES 


2x200 




8600 




74079 


12.jun.2006 


New UVES 


2x400 


35 000 


3460 


85 




12.jun.2006 


New UVES 


2x150 




5800 






12.jun.2006 


New UVES 


2x150 




8600 




77946 


ll.jun.2006 


New UVES 


2x800 


35 000 


3460 


65 




ll.jun.2006 


New UVES 


2x300 




5800 






ll.jun.2006 


New UVES 


2x300 




8600 




80003 


12.jun.2006 


New UVES 


2x2400 


35 000 


3460 


30 




12.jun.2006 


New UVES 


2x900 




5800 






12.jun.2006 


New UVES 


2x900 




8600 




80837 


04.jun.2006 


New UVES 


4x300 


35 000 


3460 


110 




04.jun.2006 


New UVES 


4x100 




5800 






04.jun.2006 


New UVES 


4x100 




8600 




81170 


05.aug.2003 


Archive 


3x1333 


40 000 


3460 


45 


85963 


04.jun.2006 


New UVES 


2x400 


35 000 


3460 


70 




04.jun.2006 


New UVES 


2x150 




5800 






04.jun.2006 


New UVES 


2x150 




8600 




87693 


10.apr.2000 


Archive 


2x3600 


40 000 


3460 


70 


88010 


07.aug.2003 


Archive 


3x1333 


40 000 


3460 


70 


92781 


10.apr.2000 


Archive 


2x1800 


40 000 


3460 


100 


94449 


12.jun.2006 


New UVES 


2x1000 


35 000 


3460 


90 




12.jun.2006 


New UVES 


2x400 




5800 






12.jun.2006 


New UVES 


2x400 




8600 




98020 


ll.apr.2001 


Archive 


2700 


50000 


3460 


110 


98532 


09.apr.2000 


Archive 


2x900 


50000 


3460 


140 


100568 


12.jun.2006 


New UVES 


2x800 


35 000 


3460 


120 




12.jun.2006 


New UVES 


2x300 




5800 






12.jun.2006 


New UVES 


2x300 




8600 




100792 


25.jun.2006 


New UVES 


2x700 


35 000 


3460 


50 




25.jun.2006 


New UVES 


2x300 




5800 






25.jun.2006 


New UVES 


2x300 




8600 




101346 


10.apr.2000 


Archive 


2x1500 


40 000 


3460 


85 


103498 


04.jun.2006 


New UVES 


2x600 


35 000 


3460 


80 




04.jun.2006 


New UVES 


2x200 




5800 






04.jun.2006 


New UVES 


2x200 




8600 




104660 


30.jun.2006 


New UVES 


2x500 


35 000 


3460 


40 




30.jun.2006 


New UVES 


2x200 




5800 






30.jun.2006 


New UVES 


2x200 




8600 




105858 


07.dec.2002 


UVES POP 


3x80 


80 000 


3460 


400 


105888 


18.oct.2005 


Archive 


800 


40000 


3460 


65 


106447 


25.jun.2006 


New UVES 


2x3600 


35 000 


3460 


45 




25.jun.2006 


New UVES 


6x1100 




5800 






30.jun.2006 


New UVES 


3600 


35 000 


3460 






30.jun.2006 


New UVES 


3x1100 




8600 




107975 


30.jun.2006 


New UVES 


2x150 


35 000 


3460 


55 




30.jun.2006 


New UVES 


2x10 




5800 






30.jun.2006 


New UVES 


2x10 




8600 




108490 


12.jun.2006 


New UVES 


300 


35 000 


3460 


90 




12.jun.2006 


New UVES 


2x30 




8600 






Uy.JUl.zUUo 


In ew U V tso 




'2 < r\r\r\ 

J J uuu 


j40U 






09.jul.2006 


New UVES 


2x30 




5800 




109067 


13.oct.2005 


New UVES 


1500 


35 000 


3460 


50 




13.oct.2005 


New UVES 


2x650 




5800 




109558 


29.may.2006 


New UVES 


1800 


35 000 


3460 


55 




29.may.2006 


New UVES 


2x600 




5800 




109646 


30.jun.2006 


New UVES 


2x400 


35 000 


3460 


105 



(1) Whether the spectra are from new observations, the ESO archive data, or the UVES POP Ubrary. 



(2) The S/N per pixel of the final combined spectrum in the Be region. 
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Table 2. continued. 



T 1 1 ri 

Hlr 


Date 01 obs. 


Type 


Exp. lime 


K 




5>/N 








(s) 




(A) 


(final) 




30.jun.2006 


New UVES 


2x150 




5800 






30.jun.2006 


New UVES 


2x150 




8600 




112229 


30.jun.2006 


New UVES 


400 


35 000 


3460 


45 




30.jun.2006 


New UVES 


2x150 




8600 




114271 


17.oct.2005 


New UVES 


1200 


35 000 


3460 


120 




17.oct.2005 


New UVES 


2x500 




8600 




114962 


29.may.2006 


New UVES 


2x600 


35 000 


3460 


100 




29.may.2006 


New UVES 


2x200 




5800 






29.may.2006 


New UVES 


2x200 




8600 




115167 


ll.aug.2003 


Archive 


3x1066 


40 000 


3460 


65 


117041 


13.oct.2005 


New UVES 


1800 


35 000 


3460 


35 




13.oct.2005 


New UVES 


2x800 




5800 




G05-19 


15.oct.2005 


New UVES 


2700 


35 000 


3460 


55 




15.oct.2005 


New UVES 


3000 


35 000 


3460 






15.oct.2005 


New UVES 


2x1300 




5800 






15.oct.2005 


New UVES 


2x1400 




8600 




G05-40 


08.oct.2001 


Archive 


2x4800 


35 000 


3460 


90 




22.nov.2001 


Archive 


2x4800 


35 000 


3460 




G66-51 


02.jun.2006 


New UVES 


2x2400 


35 000 


3460 


30 




02.jun.2006 


New UVES 


2x900 




5800 






02.jun.2006 


New UVES 


2x900 




8600 




G166-37 


02.jun.2006 


New UVES 


3600 


35 000 


3460 


25 




02.jun.2006 


New UVES 


3x1100 




5800 






12.jun.2006 


New UVES 


2x3600 


35 000 


3460 






12.jun.2006 


New UVES 


3x1100 




5800 






12.jun.2006 


New UVES 


3x1100 




8600 




Lrl /U-Zl 


Ul.jun.ZUUO 


XT-,,,, TT\/H7C 

INew U VJio 


JOUU 


J J uuu 








01.jun.2006 


New UVES 


3x1100 




5800 






02.jun.2006 


New UVES 


3600 


35 000 


3460 






02.jun.2006 


New UVES 


3x1100 




5800 






02.jun.2006 


New UVES 


2x1100 




8600 






12.jun.2006 


New UVES 


3600 


35 000 


3460 






12.jun.2006 


New UVES 


3x1100 




8600 





(1) Whether the spectra are from new observations, the ESO archive data, or the UVES POP Ubrary. 



(2) The S/N per pixel of the final combined spectrum in the Be region. 
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Table 3. Atmospheric parameters, as adopted from the hterature, beryllium abundances as derived in this work, and lithium abundances, adopted 
from the literature when available or derived in this work when not. Details on the parameters and the abundances are given in the text. The 
reference papers of the atmospheric parameters and lithium abundances are also listed. 



Star 


Teff 


logg 


logg 




[Fe/H] 


Ref. 


log(Be/H) 


log(Be/H) A(Li) Ref. 




Kelvin 


Liter. 


Paral. 


km s ' 






3131 


average 



5275 
5950 
5909 
5425 
5640 
6100 
5700 
5075 
5825 
5360 
5800 
5886 
5150 
5378 
5200 
5825 
5700 
5725 
5750 
5300 
6097 
5725 
5933 
5800 
5830 
5672 
5850 
5822 
5200 
6025 
5900 
5750 
6021 
6375 
5600 
5700 
5900 
6450 
5875 
5946 
5831 
6046 
6200 
6200 
5600 
5865 
5720 
5800 
5500 
5396 
5686 
5450 
5480 
5875 

5575 

(1) Takeda & Kawanomoto ( i2005l l. (2) 
(5) Boesgaard et al. (2005 ), (6) Fulbrig 
(9) Romano et al. (,1999i l, (10) Spite et 



HIP 171 
HIP 3026 
HIP 7459 
HIP 10140 
HIP 10449 
HIP 11952 
HIP 13366 
HIP 14086 
HIP 14594 
HIP 17001 
HIP 17147 
HIP 18802 
HIP 19007 
HIP 19814 
HIP 21609 
HIP 22632 
HIP 24030 
HIP 24316 
HIP 31188 
HIP 31639 
HIP 33221 
HIP 33582 
HIP 34285 
HIP 36491 
HIP 36640 
HIP 36818 
HIP 36849 
HIP 37853 
HIP 38625 
HIP 42592 
HIP 44075 
HIP 44124 
HIP 45554 
HIP 48152 
HIP 50139 
HIP 52771 
HIP 53070 
HIP 55022 
HIP 57265 
HIP 58145 
HIP 58962 
HIP 59490 
HIP 59750 
HIP 60632 
HIP 62882 
HIP 63559 
HIP 63918 
HIP 64426 
HIP 66665 
HIP 67655 
HIP 67863 
HIP 70681 
HIP 71458 
HIP 72461 
HIP 74067 



4.10 
3.90 
4.46 
4.10 
4.40 
4.20 
4.20 
3.60 
4.20 
3.00 
4.30 
4.33 
4.50 
4.43 
3.80 
4.30 
4.20 
4.40 
4.10 
4.30 
4.09 
4.30 
4.26 
4.40 
3.64 
4.57 
4.10 
4.42 
4.40 
4.10 
4.20 
3.60 
4.44 
4.10 
4.30 
4.50 
4.20 
4.20 
4.00 
4.41 
4.36 
4.46 
4.40 
4.40 
3.70 
4.41 
4.14 
4.10 
3.80 
4.38 
4.40 
4.50 
3.10 
4.10 
4.30 



4.30 
4.10 
4.59 
4.25 
4.38 
4.25 
4.09 
3.48 
4.42 

4.29 
4.22 
4.74 

4.55 
4.41 
4.26 
4.45 
4.28 
4.55 
4.00 
4.32 
4.30 
4.39 

4.88 
4.15 
4.18 
4.57 
4.02 
4.11 

3.96 
4.06 
4.30 
4.48 
4.27 
4.05 
4.13 
3.96 

4.46 
4.26 
4.45 

4.16 
3.91 
4.10 
3.77 
4.62 
4.42 
4.55 

4.31 
4.41 



1.05 
1.40 
1.23 
0.85 
1.00 
0.95 
0.95 
1.10 
1.10 
1.20 
1.10 
1.38 
1.20 
0.84 
1.55 
1.35 
1.11 
1.30 
1.65 
0.60 
1.86 
1.25 
1.50 
1.10 
0.70 
0.90 
1.10 
1.21 
0.30 
1.20 
1.25 
0.10 
1.34 
2.25 
0.35 
1.30 
1.45 
1.80 
1.50 
1.32 
1.30 
1.34 
1.10 
1.35 
0.04 
1.26 
1.49 
1.25 
1.05 
0.92 
1.13 
0.80 
1.98 
0.40 
1.10 



-0.90 
-1.20 
-1.15 
-1.00 
-0.80 
-1.60 
-0.70 
-0.60 
-2.00 
-2.35 
-0.80 
-0.85 
-0.50 
-0.69 
-1.60 
-1.40 
-1.20 
-1.50 
-0.70 
-0.50 
-1.30 
-0.50 
-0.90 
-0.80 
-0.70 
-0.83 
-0.70 
-0.78 
-0.70 
-2.00 
-0.70 
-1.80 
-0.75 
-2.00 
-0.70 
-1.80 
-1.40 
-0.80 
-1.00 
-1.04 
-0.80 
-1.26 
-0.60 
-1.50 
-1.10 
-0.93 
-0.65 
-0.70 
-0.80 
-0.93 
-0.70 
-1.10 
-2.10 
-2.30 
-0.80 



FOO 

FOO 
NS97 

FOO 

FOO 

FOO 

FOO 

FOO 

FOO 
GS88 

FOO 
NS97 

FOO 
NS97 

FOO 

FOO 

PrOO 

FOO 

FOO 

FOO 
NS97 

FOO 
NS97 

FOO 
GS91 
NS97 

FOO 
Ed93 

FOO 

FOO 

FOO 

FOO 
NS97 

FOO 

FOO 

FOO 

FOO 

FOO 

FOO 
NS97 
NS97 
NS97 

FOO 

FOO 

FOO 
NS97 
NS97 

FOO 

FOO 
NS97 
NS97 

FOO 
McW95 

FOO 

FOO 



-11.55 
-12.13 
-11.90 
-11.68 
-11.31 
-12.28 
-11.07 

-12.63 
<-13.83 
-11.20 
-11.41 



-12.17 
-11.55 
-11.98 
-11.46 
-11.27 
-11.67 
-11.09 
-11.85 
-11.40 
-11.44 

-11.35 
-11.16 
-11.42 
-12.58 
-11.12 
-12.92 
-11.45 
-12.67 
-11.28 

-11.80 

-11.90 
-11.31 
-11.68 
-10.60 

-12.10 
-11.20 

-10.95 
-11.31 
-11.61 
-11.69 
-10.95 
-11.12 
<-13.90 

-11.36 



-11.51 
-12.11 
-11.88 
-11.64 
-11.31 
-12.26 
-11.04 

-12.58 

-11.20 
-11.41 

<-11.64 

-12.14 
-11.53 
-11.96 
-11.46 
-11.15 
-11.67 
-11.09 
-11.85 
-11.40 
-11.42 
<-12.58 
-11.33 
-11.19 

-12.58 
-11.12 
-12.86 
-11.45 
-12.67 
-11.23 
-12.55 
-11.80 

<-12.75 
-11.90 
-11.31 
-11.64 
-10.58 

<-12.50 
-12.10 
-11.20 

-10.95 
-11.31 
-11.63 
-11.50 
-10.95 
-11.12 



-11.36 



<0.80 
2.42 
2.12 
1.54 
1.00 
2.16 
1.28 

2.18 
1.17 
1.45 
1.89 

0.98 
1.43 
2.17 
2.08 
2.00 
1.95 

<0.50 
2.32 
1.10 
2.22 
1.70 
2.36 

<0.70 
1.96 
1.37 

<0.16 
2.12 
2.06 
1.96 
2.09 
2.25 

<0.40 
2.16 
2.21 

<1.31 
2.55 
1.95 
1.95 
2.51 

<1.10 
2.21 
2.19 
2.14 
1.99 
2.01 
0.90 

<1.25 

<1.34 
1.54 
1.29 
2.22 
0.95 



Charbonnel & Primas (120 05]l, (3 ) Chen et al. (12001JI, (4) T his wo rkT 
ht (2000), (7) Favata et al. ( fTW6]l, (8 ) Ryan & Deliyannis (fT995] l, 
al. (.1994J . (11) Gratton et al. (|20U0] |. 



(1) 
(2) 
(3) 
(2) 
(4) 
(2) 
(3) 

(2) 
(2) 
(3) 
(2) 

(5) 
(2) 
(2) 
(3) 
(2) 
(4) 
(4) 
(3) 
(4) 
(3) 
(6) 
(3) 
(3) 
(3) 
(3) 
(7) 
(2) 
(3) 
(4) 
(3) 
(2) 
(4) 
(2) 
(2) 
(2) 
(5) 
(3) 
(3) 
(3) 
(3) 
(2) 
(8) 
(2) 
(3) 
(3) 
(4) 
(2) 
(2) 
(3) 
(2) 
(2) 
(4) 



Table 3. continued. 
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Star 


Teff 

Kelvin 


logg 
Liter. 


logg 
Paral. 


km s ' 


[Fe/H] 


Ref. 


log(Be/H) 
3131 


log(Be/H) 
average 


A(Li) 


Ref. 


HIP 74079 


5825 


4.00 


3.92 


1.30 


-0.70 


FOO 


-11.28 


-11.28 


2.24 


(2) 


HIP 77946 


6550 


3.60 


3.70 


1.65 


-0.90 


FOO 


<-12.95 


<-12.85 


<1.45 


(6) 


HIP 80003 


5486 


4.80 


4.80 


0.00 


-0.87 


SB02 


-11.62 


-11.62 


1.25 


(5) 


HIP 80837 


5800 


4.10 


4.07 


1.15 


-0.70 


FOO 


-11.22 


-11.20 


1.53 


(9) 


HIP 81170 


5175 


4.70 


4.64 


0.30 


-1.10 


FOO 


-11.53 


-11.42 


<0.30 


(10) 


HIP 85963 


6227 


3.94 


3.75 


2.16 


-0.71 


Ed93 


-11.78 


-11.78 


<0.70 


(4) 


HIP 87693 


6175 


4.00 


3.99 


0.90 


-2.00 


FOO 


-12.77 


-12.77 


2.20 


(2) 


HIP 88010 


5200 


4.00 


4.15 


0.70 


-1.40 


FOO 


-12.40 


-12.34 


1.67 


(2) 


HIP 92781 


5650 


4.20 


4.12 


0.95 


-0.60 


FOO 


-11.14 


-11.14 


- 


- 


HIP 94449 


5625 


3.70 


3.74 


1.15 


-1.20 


FOO 


-11.50 


-11.50 


1.81 


(6) 


HIP 98020 


5325 


4.60 


4.54 


1.10 


-1.60 


FOO 


- 


- 


1.61 


(2) 


HIP 98532 


5550 


3.60 


- 


1.30 


-1.10 


FOO 


-11.43 


-11.43 


2.15 


(3) 


HIP 100568 


5650 


4.40 


4.51 


1.10 


-1.00 


FOO 


-11.98 


-11.92 


1.92 


(2) 


HIP 100792 


5875 


4.20 


4.26 


1.40 


-1.10 


FOO 


-11.97 


-11.93 


2.02 


(2) 


HIP 101346 


6000 


3.90 


3.74 


1.40 


-0.50 


FOO 


-11.52 


-11.52 


2.15 


(9) 


HIP 103498 


5894 


4.38 


4.29 


1.32 


-1.03 


Ed93 


-11.36 


-11.39 


2.00 


(11) 


HIP 104660 


5500 


3.90 


4.01 


1.15 


-0.80 


FOO 


-11.41 


-11.44 


1.10 


(2) 


HIP 105858 


6139 


4.34 


4.26 


1.57 


-0.67 


Ed93 


-11.32 


-11.32 


2.39 


(3) 


HIP 105888 


5700 


4.00 


3.99 


1.00 


-0.60 


FOO 


-11.13 


-11.13 


1.37 


(6) 


HIP 106447 


6089 


4.04 


- 


1.50 


-2.48 


SB02 


-12.88 


- 


2.23 


(5) 


HIP 107975 


6275 


3.90 


3.87 


1.50 


-0.50 


FOO 


-12.25 


-12.25 


<1.10 


(9) 


HIP 108490 


6009 


4.41 


4.30 


1.37 


-0.72 


Ed93 


-11.34 


-11.30 


2.31 


(3) 


HIP 109067 


5300 


4.30 


4.68 


0.85 


-0.80 


FOO 


-11.39 


-11.37 


<0.10 


(4) 


HIP 109558 


6025 


4.00 


4.09 


1.10 


-1.50 


FOO 


-12.30 


-12.34 


2.10 


(2) 


HIP 109646 


5910 


4.25 


4.21 


1.50 


-0.64 


Ed93 


-11.38 


-11.36 


2.23 


(3) 


HIP 112229 


5983 


4.37 


4.16 


1.40 


-0.65 


Ed93 


-11.21 


-11.25 


2.31 


(3) 


HIP 114271 


6200 


4.10 


4.21 


0.90 


-1.70 


FOO 


-12.44 


-12.39 


2.32 


(2) 


HIP 1 14962 


5825 


4.30 


3.81 


1.40 


-1.40 


FOO 


-12.40 


-12.34 


2.25 


(2) 


HIP 115167 


6100 


3.80 


3.57 


1.15 


-1.50 


FOO 


-12.45 


-12.45 


2.17 


(2) 


HIP 117041 


5300 


4.20 


4.11 


0.90 


-0.80 


FOO 


-11.34 


-11.34 


<0.10 


(4) 


G05-19 


5942 


4.24 




0.89 


-1.10 


SB02 


-11.98 


-12.01 


2.26 


(5) 


G05-40 


5863 


4.24 




1.48 


-0.83 


NS97 


-11.15 


-11.15 


1.90 


(3) 


G66-51 


5255 


4.48 




0.90 


-1.00 


PrOO 


-11.47 


-11.37 


<0.30 


(4) 


G166-37 


5300 


4.80 


4.69 


0.60 


-1.20 


FOO 






1.28 


(5) 


G170-21 


5664 


4.65 




0.00 


-1.45 


SB02 


-11.65 


-11.58 


1.91 


(5) 



(1) Takeda & Kawa nomoto (12005 ), (2) Charbonnel & Primas (2005), (3) Chen et al. (1200 iD , (4) T his wo rk, 
(5) Boesgaard et a l. (|2005] l, (6) Fulbright ( 2000) (7) Favata et al. (.1996), (8 ) Ryan & DeUyannis (fT995] l, 
(9) Romano et al. (fT999] l, (10) Spite et al. ( fT994] l, (11) Gratton et al. (|20D0| |. 
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Table 6. The probabihties of each star of belonging to the thin disk, the thick disk, and the halo as calculated by Venn i2004b . the kinematic 
classification' accor ding to Gratton et al. j2003a||2003bl l, the average [a/Fe], and the orbital parameters: minimum radii (Rrmin)i maximum radii 
(Rmax)i maximum distance from the plane (Z„ax)> eccentricity, and rotational velocity. 



Star 


1 nin uisK 


1 nicK uisK 


Halo 


Gratton 


[ff/rej 


R min. 


R max. 


Z max 


ecc. 


Vrot 




prob. 


prob. 


prob. 


clas. 














HIP 171 


0.20 


0.80 


0.00 





0.41 


4.58 


8.50 


0.31 


0.30 


152.63 


HIP 3026 


0.00 


0.00 


1.00 


1 


0.30 


0.19 


10.72 


7.62 


0.97 


-10.61 


HIP 7459 


0.00 


0.00 


1.00 


1 


0.17 


0.67 


16.87 


8.34 


0.92 


-38.62 


HIP 10140 


0.00 


0.90 


0.10 





0.29 


4.29 


8.97 


0.61 


0.35 


147.98 


T TTTl -I f\ A A r\ 

HIP 10449 


0.00 


0.00 


1.00 


1 


0.33 


0.40 


12.82 


6.64 


0.94 


22.69 


HIP 11952^ 


0.00 


0.90 


0.10 





0.43 


3.70 


8.67 


0.42 


0.40 


122.0 


HIP 13366 


0.00 


0.70 


0.30 





0.31 


3.35 


8.85 


1.65 


0.45 


120.06 


HIP 14086 


0.10 


0.90 


0.00 





0.34 


5.27 


9.02 


1.05 


0.26 


167.76 


T TTT^ 1 A C r\ A 

HIP 14594 


0.00 


0.10 


0.90 





0.41 


2.34 


12.24 


1.52 


0.68 


102.15 


HIP 17001 


0.00 


0.00 


1.00 




0.16 










— 


HIP 17147 


0.00 


0.90 


0.10 


' i 


0.33 


V ■ ■ 


9.93 


0.52 


0.45 


140.82 


HIP 18802 


0.30 


0.70 


0.00 





0.18 


4.30 


11.06 


0.05 


0.44 


162.33 


HIP 19007 


0.90 


0.10 


0.00 


2 


0.21 


7.95 


10.97 


0.11 


0.16 


236.72 


HIP 198 14^ 


0.00 


0.00 


1.00 


1 


0.05 


0.58 


31.41 


10.39 


0.96 


28.00 


HIP 21609 


0.00 


0.00 


1.00 


1 


0.45 


1.51 


47.65 


2.24 


0.94 


89.89 


HIP 22632 


0.00 


0.90 


0.10 





0.36 


3.48 


8.87 


0.24 


0.44 


128.43 


HIP 24030 


0.00 


0.80 


0.20 





0.28 


4.06 


8.66 


2.49 


0.36 


131.64 


HIP 24316 


0.00 


0.00 


1.00 


1 


0.37 


2.65 


15.78 


4.71 


0.71 


-120.15 


HIP 31188 


0.00 


0.90 


0.10 





0.13 


6.98 


8.68 


1.42 


0.11 


191.26 


HIP 31639 


0.20 


0.70 


0.00 





0.27 


5.71 


8.55 


0.55 


0.20 


173.93 


HIP 33221 


0.00 


0.70 


0.30 





0.28 


6.67 


8.88 


3.43 


0.14 


175.82 


HIP 33582 


0.00 


0.00 


1.00 





0.37 


1.74 


14.68 


0.18 


0.79 


89.37 


HIP 34285 


0.00 


0.00 


1.00 


1 


0.05 


0.63 


25.25 


2.14 


0.95 


-43.26 


HIP 36491 


0.00 


0.80 


0.20 





0.30 


2.60 


8.71 


0.04 


0.54 


103.85 
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0.00 


0.00 


1.00 




0.06 












HIP 36818 


0.00 


0.00 
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1 


0.04 


0.66 


14.63 


1.73 


0.91 


-37.34 
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0.00 


0.90 


0.10 





0.26 


3.69 


9.27 


0.78 


0.43 
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0.00 


0.00 
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0.28 


3.96 


11.00 


0.38 


0.47 
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0.20 


0.80 


0.00 





0.28 


4.92 
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0.08 


0.31 


166.78 
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HIP 57265 


0.00 


0.00 


1.00 


1 


0.21 


0.13 


43.31 


21.72 


0.99 


-6.00 


HIP 58145 


0.00 


0.90 


0.10 





0.24 


4.87 


8.55 


1.37 


0.27 


143.00 


HIP 58962^ 


0.00 


0.00 


1.00 


1 


0.06 


0.56 


8.54 


2.75 


0.88 


21.00 


HIP 59490 


0.00 


0.00 


1.00 


1 


0.25 


0.82 


8.61 


0.55 


0.83 


-44.47 


HIP 59750 


0.00 


0.90 


0.10 





0.31 


4.36 


9.23 


0.89 


0.36 


150.82 


HIP 60632 


0.00 


0.00 


1.00 


1 


0.41 


0.21 


10.92 


6.77 


0.96 


11.03 


HIP 62882 


0.00 


0.00 


1.00 


1 


0.37 


0.33 


21.76 


10.23 


0.97 


22.14 


HIP 63559 


0.00 


0.00 


1.00 


1 


0.11 


0.38 


9.00 


5.19 


0.92 


-18.24 
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HIP 64426 


0.00 


0.90 


0.10 





0.29 


4.58 


9.52 


1.35 


0.35 


155.72 


HIP 66665^ 


0.00 


0.20 


0.80 





0.36 


3.28 


12.57 


2.19 


0.59 


123.00 


HIP 67655 


0.30 


0.70 


0.00 





0.22 


5.82 


8.75 


0.38 


0.19 


177.72 


HIP 67863 


0.00 


0.00 


1.00 


1 


0.22 


0.59 


9.01 


6.24 


0.88 


-27.24 


HIP 70681 


0.00 


0.90 


0.10 





0.31 


6.14 


8.50 


1.35 


0.16 


178.04 


HIP 71458 


0.00 


0.00 


1.00 




0.35 












HIP 72461 


0.00 


0.60 


0.40 





0.44 


2.83 


10.82 


0.46 


0.59 


119.11 


HIP 74067 


0.00 


0.90 


0.10 





0.30 


4.95 


8.53 


1.09 


0.27 


157.5 



(1) The following code is adopted: for the accretion component , 1 for the dissipative component, and 2 for thin disk stars. 



(2) The orbital parameters for these stars were calculated in this work. 
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Table 6. continued. 
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0.00 
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1 
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3.19 
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4.32 


0.78 
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0.00 


0.00 
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1 


0.32 


0.70 
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0.87 
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0.00 
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0.36 
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0.72 
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0.00 


2 
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0.05 


0.11 
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0.00 


1.00 


1 


0.44 
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(1) The following code is adopted: for the accretion component , 1 for the dissipative component, and 2 for thin disk stars. 



(2) The orbital parameters for these stars were calculated in this work. 



Appendix A: Detailed comparison 

In this section we present a detailed comparison of our Be abundance results with previous results from the literature on a star-by- 
star basis. In most of the cases where differences in the abundances are seen, differences in the log g values can also be found. In 
these cases we are confident our results are robust since our log g values show excellent agreement with the ones derived using 
Hipparcos parallaxes. 

HIP 171 (HD 224930) 

The Be abundance of star HIP 171 was previously determined by Stephens et al. d 19971 1. Iog(Be/H) = -1 1.02. The value we found, 
Iog(Be/H) — -1 1.55, is somewhat different. This difference is most probably a result of the different log g adopted; log g - 4.62 by 
Stephens et al. (1 19971 1 and log g = 4. 10 by us. 

HIP 11952 (HD 16031) 

The Be abundance of HIP 1 1952 was previously determined by Gibnore et al. d 19921 ). Adopting log g - 3.90 and [Fe/H] = -1.96 
they found Iog(Be/H) - -12.37. Adopting log g = 4.20 and [Fe/H] - -1.60, we found log(Be/H) = -12.28, which shows a good 
agreement within the uncertainties, in spite of the different atmospheric parameters. 
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HIP 14594 (HD 19445) 

The Be abundance of star HIP 14594 was previously determined by four papers. Rebolo et al. d 19881 1 found an upper limit of 
log(Be/H) < -11.70, Ryan et al. (1990") found an upper limit of log(Be/H) < -12.30, Boesgaard & King (^1993 ) found log(Be/H) = 
-12.14, and Boesgaard et al. (1999 ) found a range of values from log(Be/H)= -12.45 to -12.55. The value we derived, log(Be/H) 
= -12.63, is in good agreement with the lower limit of the range derived by Boesgaard et al. ( I1999I ). 

HIP 17001 (CD-24 1782) 

The Be abundance of star HIP 17001 was previously determined by Garcia Perez & Primas ( I2006I I who found log(Be/H) - -13.45 
in LTE and log(Be/H) = -13.54 in NLTE. We derived an upper limit of log(Be/H) < -13.83. The difference in the values is mainly 
due to the different log g adopted, 3.00 by us and 3.46 by Garcia Perez & Primas ( 12006! ). 

HIP 17147 (HD 22879) 

The Be abundance of HIP 17147 was previously determined by Beckman et al. ( 119891 1. log(Be/H) = -1 1.25, This value agrees well 
with the one found in this work, log(Be/H) = -11 .20. 

HIP 18802 (HD 25704) 

The Be abundance of star HIP 1 8802 was previously determined by Molaro et al. |1997a[ log(Be/H) = - 1 1 .6 1 . In this work we found 
log(Be/H) = -11.41. The log g values, although similar, log g - 4.20 by Molaro et al. |1997al and log g = 4.33 by us, account for a 
difference of ~ 0.08 dex in the abundance, what would bring them to an agreement within the uncertainties. 

HIP 24316 (HD 34328) 

The Be abundance of HIP 24316 was previously determined by Gilmore et al. ( 119921 ), log(Be/H) - -1 1.90. This value agrees well 
with the value found in this work, log(Be/H) = -1 1.98. 

HIP 37853 (HR3018) 

The Be abundance of HIP 37853 was previously determined by Gilmore et al. ( 119921 ). log(Be/H) - -1 1.20. This value agrees well 
with the one found in this work, log(Be/H) = -11.16. 

HIP 42592 (HD 74000) 

The Be abundance of star HIP 42592 was previously determined by Ryan et al. (119901) . log(Be/H) < -12.20 and by Boesgaard et 
al. ( 119991 1. which found a range of values from log(Be/H) - -12.03 to -12.47. The value found in this work, log(Be/H) - -12.58, 
agrees with the lower limit of the range determined by Boesgaard et al. ('1999'), within the uncertainties. The small difference can be 
removed if we adopt the same log g used by Boesgaard et al., 0. 16 dex higher than ours (increasing our abundance by ~ 0.08 dex). 

HIP 44075 (HD 76932) 

The Be abundance of star HIP 44075 was determined by a number of papers in literature. Molaro & Beckman ( 119841 ) found an upper 
limit of log(Be/H) < -1 1.52, Beckman et al. (1989) fo und log(Be/H) = -1 1.69, both Gilmore et al. ( [1992]) and Ryan et al. ( 199% 
found log(Be/H) = -11.30, Boesgaard & Kin g ( [T993]) found log(Be/H) = -11.04, Gar cia Lop ez et al. ( 11995bl ) found log(Be/H) = 
-11.36 (in NLTE), Thorburn & Hobbs ( fT996l ) found log(Be/H) = -1 1.45, Molaro et al. ( |1997al ) found -1 1.21, and Boesgaard et al. 
(119991 ) found a range of values from log(Be/H)= -11.17 to -11. 24. The value we found, log(Be/H) = - 1 1 . 1 2, is in the upper range 
of the values listed above and in good agreement with Boesgaard et al. (119991 ). 

HIP 48152 (HD 84937) 

For star H IP 481 52, Ryan et al. ( [1992] ) found log(Be/H) < -12.85 , Boes gaard & King ([1993] ) found log(Be/H) = -12.85, Thorburn 
& Hobbs ([19961) found log(Be/H) < -12.95, and Boesgaard et al. ([T9W ), from log(Be/H)= -12.83 to -1 2.94. The value we derived 
is log(Be/H) - -12.67. Reducing our log g by 0.20 dex to match the one adopted by Boesgaard et al. (119991 ). for example, would 
reduce our abundance by ~ 0.1 1 dex, resulting in an agreement within the uncertainties. 

HIP 53070 (HD 94028) 

For HIP 53070, Boesgaard & King (T993^ found log(Be/H) = -11.56, Garcia Lopez et al. (T995F) found log(Be/H) = -11.66 
(in NLTE), Thorburn & Hobbs (1996) found log(Be/H) = -11.65, and Boesgaard et al. (1999 ) found a range of values from 
log(Be/H)= -1 1.51 to -1 1.55. The value we derived is log(Be/H) - -1 1.80. Once again, a change in the adopted log g would bring 
the abundances to a better agreement. 
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HIP 55022 (HD 97916) 

For HIP 55022, Boesgaard (l2007b found log(Be/H) < -13.30 while we derived log(Be/H) < = -12.75. 
HIP 59750 (HD 106516) 

For HIP 59750, Molaro et al. ( !1997al l found log(Be/H) < -12.76 and Stephens et al. (T997) found log(Be/H) < -12.61. The upper- 
limit we derived is log(Be/H) < -12.50, agrees with the previous results within the uncertainties. 

HIP 64426 (HD 114762) 

For H IP 644 26, Stephens et al. ( fT997 ^ found log(Be/H) = -11.05, Boesgaard & King (|T993T l found log(Be/H) = -11.14, Santos 
et al. ( [2002] l found log(Be/H) = -11.03, and Santos et al. (2004) found log(Be/H) = -11.18. The value we derived is log(Be/H) = 
-11.31. The difference is again related to the different choice of log g values. 

HIP 71458 (HD 128279) 

For HIP 71458, Molai'o et al. ( !1997al l found log(Be/H) = -12.75 and by Garcia Perez & Primas ( |2006l ) found log(Be/H) < -14.01, 
in LTE, and log(Be/H) = -13.94 in NLTE. We derived an upper limit of log(Be/H) < -13.90. 

HIP 74079 (HD 134169) 

The Be abundance of HIP 74079 was previously determined by five papers. Gilmore et al. (|1992| l and Ryan et al. (1X992') both found 
log(Be/H) = -1 1.35, Boesgaard & King a993) found log(Be/H) = -1 1.29, Garcia Lopez et al. ( |1995bl ) found log(Be/H) = -1 1.23 
(in NLTE), and Boesgaard et al. ( 199?ll found a range of values from log(Be/H) — -11.32 to log(Be/H) — —11.40. All these values 
agree with ours, log(Be/H) = -1 1.28, within the uncertainties. 

HIP 80837 (HD 148816) 

The Be abundance of HIP 80837 was previously determined by Stephens et al. (119971 ), log(Be/H) = -11.07 and by Boesgaard &. 
King ( I1993I I, log(Be/H) - -11.10. Both values are consistent with ours, log(Be/H) - -1 1.22, within the uncertainties. 

HIP 87693 (BD+20 3603) 

The Be abundance of star HIP 87693 was previously determined by Boesgaard et al. ( 11999! ) who found a range of values from 
log(Be/H) - -12.40 to log(Be/H) - -12.62. Our value, log(Be/H) - -12.77, is again consistent, within the uncertainties, with the 
lower limit of the range found by Boesgaard et al. ( 119991 1. The difference in the adopted log g values, 4.33 or 4.27 by Boesgaard et 
al. and 4.00 by us, seems to be the main reason for the different abundances. 

HIP 98532 (HD 189558) 

The Be abundance of star HIP 98532 was previously determined by Rebolo et al. ( 11988! ) who found log(Be/H) = -11.70 and by 
Boesgaard & King ("1993) who found log(Be/H) - -10.99. Both papers claim higher uncertainties for this star when compared 
with the other stars of the sample. We determined a value of log(Be/H) = -1 1.43, which is intermediate between the two previous 
determinations. 

HIP 100792 (HD 194598) 

The Be abundance of star HIP 100792 was previously determined by three papers. Rebolo et al. ( 11988! ) found log(Be/H) = -1 1.70, 
Thorburn & Hobbs ( fT996t found log(Be/H) = - 1 1 .95, and Boesgaard et al. ( fT999l ) found a range of values from log(Be/H)= -11.73 
to -11.88. The value determined in this work is log(Be/H) = -11.97. Our result agrees very well with the one from Thorburn & 
Hobbs ( 1996) although they adopt a smaller log g, 4.00, when compared to ours, log g - 4.20. Within the uncertainties, our result 
agrees with the lower limit of the range of values found by Boesgaard et al. ( |1999! l. 

HIP 101346 (HD 195633) 

The Be abundance of star HIP 101346 was previously determined in three papers. Boesgaard & King (119931) found log(Be/H) = 
-1 1.21, Stephens et al. ( [19971 ) found log(Be/H) = -1 1.29, and Boesgaard & Novicki (2006 ) found log(Be/H)= -1 1.34. Our value is 
somewhat smaller than these, log(Be/H) = -11 .52. The parameter controlling most of this difference seems to be the metallicity. We 
adopt [Fe/ H] = -0 .50, while Boesgaard & King ( [1993] ) adopt [Fe/H] = -1.07, Stephens et al. ( 1997) [Fe/H] = -1.00, and Boesgaard 
& Novicki (12006! ) [Fe/H] - -0.88. Adopting the value of Boesgaard & Novicki ( l2006! l. for example, would increase our result by 
0.08 dex, bringing the values to agree within the uncertainties. 
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HIP 104660 (HD 201889) 

The Be abundance of star HIP 104660 was previously determined by Boesgaard & King ( |1993l l. log(Be/H) - -11.43, and by 
Boesgaard et al. (1999') which found a range of values from log(Be/H)= -11.30 to -11.38. Our value, log(Be/H) = -11.41, is in 
excellent agreement with these. 

HIP 105858 (HR8181) 

The Be abundance of HIP 105858 was previously determined by Gilmore et al. ( 119921 ). A value of log(Be/H) - -1 1.40 was found, 
in excellent agreement with ours, log(Be/H) - -1 1.32. 

HIP 107975 (HD 207978) 

The Be abundance of star HIP 1 07975 was previously determined by Stephens et al. ( 119971 1. which found an upper limit of log(Be/H) 
< -12.38. Our analysis provides a detection mainly based on line 3131 A, log(Be/H) = -12.25. 

HIP 108490 (HD 208906) 

The Be abundance of star HIP 108490 was previously determined in three papers. Boesgaard & King ( I1993I I found log(Be/H) = 
-11.10, Stephens et al. (TWT) found log(Be/H) = -11.19, and Boesgaard et al. (200^) found log(Be/H)= -1 1.30. This last value is 
in excellent agreement with our own, log(Be/H) = -1 1.34. The difference with Boesgaard & King (1993 ) can be explained by the 
smaller log g adopted by them, log g = 4.25, when compared with the one adopted in this work, log g - 4.41 (and also by Boesgaard 
et al. 2004). 

HIP 109558 (BD+17 4708) 

The Be abundance of star HIP 109558 was previously determined by Boesgaard et al. ( 119991 ), who found a range of values from 
log(Be/H)= -12.28 to -12.42. Our result, log(Be/H) - -12.30, is in excellent agreement with this range. 

HIP 114271 (HD 218502) 

The Be abundance of HIP 1 14271 was previously determined by Molaro et al. ( |1997al l, log(Be/H)= -12.56. Our result, log(Be/H) 
= -12.44, agrees with it within the uncertainties. 

HIP 114962 (HD 219617) 

For H IP 1 14 962, Rebolo et al. ( fT988] l found log(Be/H) < -1 1.60, Molaro et al. ( !1997aD found log(Be/H) = -12.56, and Boesgaard 
et al. ( 11999! ) found a range from log(Be/H) = -12.09 to -12.15. Our value, log(Be/H) - -12.40, is closer to the value found by 
Molaro et al. ( |1997al l than to the one by Boesgaard et al. d 19991 1 although our parameters are very different from the former and 
similar to the latter. 

Appendix B: FOO subsample 

In this appendix we show plots of the relations between log(Be/H) and [Fe/H] (Fig. IB. II ) and between log(Be/H) and [c/H] (Fig. 
|B.2| i for the stars of the FOO subsample. The fits are statistically identical to the ones obtained with the whole clean sample. 

Appendix C: Accretion and Dissipative Components 

In this appendix we divide our sample in an "accretion" and a "dissipative" component according to the prescriptions of Gratton et 
al. (I2003all2003bl . 

The purpose is to test whether the same characteristics of our sample would be obtained with this different classification. The 
division of the stars according to this classification is also listed in Tabled Some of the stars classified as thin disk stars by Venn et 
al. C2004j are classified as stars of the dissipative component by the criteria of Gratton et al. ( |2003al|2003b ). In the clean sample we 
have 34 stars classified as accretion component and 39 as dissipative component. In figures ICTl IC.2I ICJl and lC.4l we show how 
these components behave in each of the plots presented before for the halo and thick disk components. Therefore, these should be 
directly compared with figures [T2l[l9l|20] As it can be seen, the overlap is really high and the conclusions obtained are the same. 
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Fig.B.l. Diagram of [Fe/H] vs. log(Be/H) for the stars analyzed by Fulbright (I2000I I. Only detections are shown. An example of 
error bar is shown in the lower right corner. 
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Fig.B.2. Diagram of [a/H] vs. log(Be/H) for the stars analyzed by Fulbright (120001 ). Only detections are shown. An example of 
error bar is shown in the lower right corner. 
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Fig. C.l. Diagram of [a/Fe] vs. log(Be/H). The dissipative component stars are shown in the upper panel while the accretion 
component stars are shown in the lower panel. The dissipative stars diagram is characterized by a large scatter while the accretion 
stars diagram clearly divides in two sequences. 
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Fig. C.l. Diagram of log(Be/H) vs. the V, the component of the space velocity of the star in the direction of the disk rotation. The 
dissipative component stars are shown as filled circles, accretion stars are shown as starred symbols, and the subgroup of accretion 
stars as open circles. A typical error of + 12 Km s"' in V was adopted (see Gratton et al. l2003b] l. 
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Fig. C.3. Diagram of log(Be/H) vs. R min, the perigalactic distance of stellar orbit. The dissipative stars are shown as filled circles, 
accretion stars are shown as starred symbols, and the subgroup of accretion stars as open circles. A typical error of + 0.40 Kpc in 
Rmin was adopted (see Gratton et al. |2003b] l. 
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Fig. C.4. Diagram of log(Be/H) vs. [Fe/H] where only the accretion stars are shown. The subgroup of stars with low alpha is shown 
as open circles, the remaining accretion stars are shown as starred symbols. The linear fit for all the accretion stars is shown to guide 
the eye. 



